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A8STHACT 


Ifcny  writers  have  been  ccncer’^ed  with  the  problaa  of  justtf/lng 
the  vme  of  ste^dard  tests  baised  on  norrnAl  theory  essua^tiona  when  io 
fact  the  underlying  distribution  Is  suspected  of  being  non-nomiJL. 
ntle  research  has  been  concerned  sainXy  with  testa  based  only  on 
vslvmrlate  situations .  They  have  found  that  the  stsundard  t  and 
f  tests  are  restarlcably  Ineenatblve  to  deviations  from  nonaal.lty. 

Here,  the  nultivarlate  jianerall cation  of  the  t  test  is  lovesti- 
4Bste4.  The  first  four  pensutstlos  cusMlaots  are  dsterniEMd  for  a 
•tatlstlc  which  is  a  slx^tle  function  of  tbs  cr—miljr  used  Hotf.ULlsc'* 
T*.  In  the  wlwmrlate  situatioo  T*  slsgiljr  heeosws  the  Siiuare  of 
“b,  and  tests  haeed  on  ?*  are  Identical  with  tests  beeed  on  t. 
Cavlete  detslls  of  the  Investlcatloo  are  carried  tbroagb  for  only 
%«o  dlJMaalaae  Irat  sixilsur  Methods  should  tpply  to  aore  diJHBSioae . 
XlH  derlTsd  peraKCatlos  eueulMte  are  applied  to  date  oMwlaad  txbm 
»  saagnUiig  aaperlMot.  ns  seaplse  are  fro«  hi  variate  Boiml  and 
(/ raetstfidwlar  povulatlone.  Rm  taaplse  drawn  from  tha  noraal  popula- 

tleo  aervad  aa  a  verlficatloo  of  the  theoretical  coagnitatloos  and 
also  providsd  a  oaaparlsoo  for  the  sapirioal  results  baaed  on  the 
I  aaapXee  froa  the  reeteumilar  population.  The  actual  esqjtirlcskL  cal- 
eulatiois  consisted  in  obtaining  all  the  perautatioos  of  the  statistic 
umtd  (a  fueetion  of  for  each  ssusple  and  obtaining  freQuea<^  dls- 

tvlstisluiM  of  the  or  pereutatloios  which  give  values  for  the 

statistic  whj  ch  are  greeter  thsu>  tt.o«-e  obtained  frews  aoraal  theory 
aaBiagstlans  for  a  certain  aigoificance  level.  It  is  observed  that 


there  la  Blld  dls&greeaent  betweeti  the  noainaLl  percentage  point  of 
the  teat  bued  oa  the  asauigitloo  that  the  underlying  dlstrlbutlcm  la 
aoiaal  and  the  actual  percentage  point  obtained  by  caoBidering  all 
perautatlona  of  each  aajqxle  froa  the  rectangular  population. 

K  dlacuaaloo  of  varloua  vaya  of  adjusting  the  teat  criterion  Is 
included  whan  coe  auapeeti  that  the  data  Coat  not  come  froa  a  noraal 
diatrlbotiod.  Thaae  aathodji  are  used  mi  the  8aiq>lea  from  the  rect¬ 
angular  dlatrlbutloo .  0n«  of  theae  intboda  proeldea  excellant  agree- 

aent  at  the  varloua  algnlfleaoce  levela  loveatl gated. 

Sogaaatlcna  toward  dlreetlona  for  future  Inveatlgatlooa  are  given. 
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Chapter  1 

IKTRODUCTIOK^  AKD  ORS-DIMBISIOKAX.  CASE 

1.1  tatroductloo 

Tba  probleE  ot  JmtiTjtetg  or  aodlf^flnc  tanalysis  of  vmrlMiet  pro- 
oedurei  vhea  the  font  of  the  populAtloB  gowmlng  the  obeerrwtlooa  le 
Dot  knowB  to  bo  Doroal  has  beoa  invaetisated  by  Pltaea  [12]>  Welch  [lb] 
and  otbera  [1,  3,  7>  d,  9}.  The  earlier  papers  of  Pitoan  [10,  U] 
aada  It  clear  that  be  orlslnally  set  out  to  find  an  i^roxlaately  noo- 
paraaetric  test  based  on  quadratic  exprsasioBs ,  Pltaas  [12]  •oes 
deeper  Into  the  tandeed  zatloe  theory  for  the  sliqpla  two-way  layout 
tfakj  does  Welch. 

nw  problea  considered  bare  is  the  extensioe  of  Htima’e  anelysls- 
of -variance  laodoalattlOQ  test  to  the  case  of  sera  than  oeo  respoase. 

Such  an  exteaeloo  aij^t^  ie  tho  case  of  two  traataente,  be  expected  to 
lead  to  Justlftcatloo  or  sodlflcatlon  of  the  taae  of  BataUlDK'e 
(9).  full  detail#  bare  bees  oarrled  tlu^Aidh  for  ceily  the  caee  of  tvo 
treataeots  and  two  respoossa. 

2.  yitaan*#  great— it  of  the  Onlwailete  Cece 

Slnca  the  ■nltlvarlste  pxtAlea  vlll  be  treated  aluoc  the  lines  of 
PltMa*e  treetaent  of  the  univariate  csee,  it  ueaos  ioflcal  to  begin  Igr 
CSitlialBg  his  pcroeedurs.  Be  eoBstdere  k  blocks  (called  '^tidies*  by 
yttaan),.  each  eoBniattes  sf  :  icdXTldusle.  mCq  »et  of  n  lodivldMals 
is  Bot  aeeesearlly  considered  as  having  cone  froas  a  larger  pepulaticR. 
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n*  D  traataents  ara  aaclgnad  tc  tha  IndlrlduaL*  ot  a  block  at  raadon. 
Tte  exparloental  arrssssssst  Is  thus  la  tiie  fom  of  randctalsafl  blocks 
and  tba  lotaraat  la  la  tasting  idtether  dlffareuues  In  treataenta  hava 
produead  aajr  rtal  dlffaraacas  In  tba  r^t^poase  measurad. 

For  k  blocks  1  aaeta  vlth  n  treataienta,  tha  laspoosas  x^j  can 
ba  dlaplajad  as  foUova 


*u 

*3.2 

•••  »ifl 

•*21 

*22 

•**  *2a 

*kl 

*k2 

•••  *ka 

*lr'  *2r'  •••'  *kr 

ara  tha> 

rasponaas  of  tba  k  ladlviduals 

atibjactad  to  traatstaat  r. 

Lst 


1. 


a  j5l  *1J  ' 

k  l5l  *1J  ' 


1.1 


Md  X  ■  ^  .-1  3t  .  “  cr  .^1  .“i  ».  • 

..  a  .J  c  14  1.  ka  1h1  J<d.  IJ 

k  B  ^ 

Tte  total  SUB  of  avtaraa  3  1<<  dafioad  to  be  j£l  ^ij.  nis 


of  agMuaa  eaa  ba  partltini^  »  foUoss: 


k  n 

ill  A  '*ij  ”  *1.  *  *.j  ■*■ 


aliera  tba  four  ^juaiititlas  tm  the  rl^t  of  aquBtica  1.2  ara  raapaetlvaly 
tha  svBi  of  e^:i£azaa  (N)  aaalgaabla  to  the  ovamll  aaar,  blocks  (a), 
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treataenta  (B)  aod  rasldu&I  (C). 

The  usuail  analyala  of  varlanca  ratio  for  testing  the  null  hypotheala 
Ro!  the  effect  asalgnahle  to  traetment  j  is  zero;  versus  the  altenui' 
tlva  hypothesis  Ht.:  the  effect  assignable  to  treatatent  J  is  not  zero/ 
for  all  J  Is 

A  _  ^  3 

A  A 


In  order  for  this  ratio  to  have  aa  r  distribution  with  (n>l)  and 
(k«l}(a>l)  dagreoa  of  freedoa,  it  la  sufficiont  that  the  be 

raodon  lodepeodent  observations  froa  normal  populations  with  the  sama 
vsolaoces.  It  is  also  necessary  to  assutae  that  the  effects  of  blocks 
and  treataants  <u*  additive.  Ana!grtically,  this  last  assuaption 
states  that  the  aeans  of  the  nornaX  population  associated  with  the  Indi¬ 
vidual  ealla  are  assuned  to  be  of  the  fora: 

“IJ  - 


with 


ZBj  o  0  and  •  0. 


n»  panuietar  p  is  the  average  of  all  the  population  means,  is 

the  average  effect  assignable  to  block  i.  The  sssusptlon  of  additi¬ 
vity  implies  that  if  block  i  gives  rise  to  a  masureaent  ten  units 
greater  than  the  measureaent  on  the  aaoe  treatment  in  some  other  block, 
then  any  other  jneasuresient  In  block  1  will  be  ten  unite  larger  (in 
the  population  mean)  liUin  the  corresponding  messurement  In  the  other 


block , 
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It  c«a  te  sbovn  that  tb«  eujbs  of  (qvarei  on  the  rlglit  hand  side  of 
•quation  1.2,  vhe»  dlvidad  by  o^,  the  cognon  population  variance, 
are  each  icdependentljr  dlatributed  by  ohl>squaa«  lavs  vith  (k'l)(n-l), 
(k-i)  aad  (n*l)  and  1  dagresR  of  freedom  respectively.  It  foUowa, 
then,  that  under  suitable  assuicptloDa  expression  1.3  is  distributed 
according  to  un  F  vltb  (n-l)  and  (K-l)(a--l)  degrees  of  freedom. 
Instead,  of  the  f  ratio,  Pitman  (12]  considered  the  ratio 

„  ECTreataenta)^ _ 

"  iXTreat^nta)*  +  (CiReaiciuai)^ 

This  ratio  will  now  be  denoted  by  W  ■  B/(B  +  C).  Under  the  above 
aasuaptlaas  B  and  C  when  divided  by  cr^  ars  lndependentl.y  distri¬ 
buted  as  X^'s  with  (n-l)  and  (k-l)(n-l)  degrees  of  freedom  respec¬ 
tively.  It  foUowB  that  W  has  a  beta  dlatrlbutioo  with  (n-l)  and 

(k-l)(n-l)  degrees  of  freedoa. 

Pltaan  [12]  treated  the  pit>bles  of  testing  the  null  hypothesis 

that  ths  trsataents  are  all  equal  without  making  any  assueptlon  about 

the  X,,.  If  the  null  hypothesis  Is  In  fact  true,  the  value  of  W 
*  J 

observed  is  ths  result  of  the  treatsegnte,  now  mens  labels,  being  dis¬ 
tributed  at  randan  to  the  various  lodlvlduals  in  the  blocks. 

An  eqptally  llkaiy  value  for  W  would  urlse  if  the  observed  res- 
poneea  were  shuffled  within  the  blocks.  Consider  for  axastple  the  case 
of  three  treatmta  and  thvM  blocks  vhnrs  the  data  ere  given  in  Table 
1.1(a).  The  valva  of  if  frost  this  set  of  obse'.va-loos  Is  21/6?. 

It  the  null  hypothesis  Is  true  thee  the  observed  vaaues  4,  2,  -3  in 
Block  1  have  appeared  in  the  order  2,  -3,  k,  eay,  and  Blocks 
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?  «jnd  3  might  have  rea-i  1,  -2,  2  and  0,  -ij  3  respectively  as 
given  in  Table  1.1(b).  After  this  ehuffilug,  W  »  5?/65.  It  Is  to 
ba  noted  that  the  block  sums  of  squares  remalne  unchanged  under  such 
reshuffling.  The  sums  of  squares  which  do  vary  are  treatak'nt  and 


Tahls  1.1 

Treotraents  Treattnents 


Blocks  1 

2 

i 

Blocks 

1 

2 

i 

1  k 

2 

-3 

1 

2 

-3 

If 

Oi 

<MI| 

-2 

1 

2 

1 

-2 

2 

1  0 

3 

-i 

i 

0 

3 

Treatiect  sum  of  squsuvs ; 

14 

38 

Besldual  sum  of  squares; 

88/3 

16/3 

(*) 

(b) 

residual.  The  nuiA)er  of  different  values  that  U  csui  take  oc  for 
this  cxa^>le  Is  exactly  the  nw^r  of  ways  In  which  the  blocks  can  be 
reshuffled  in  such  a  way  that  the  treataent  aeons  are  different  frost 
every  other  set  of  trestseut  aeaos  obtained  frost  the  other  shufflings 
Hie  nusiber  of  values  of  W,  each  value  being  equally  likely,  for  the 
case  of  n  treatments  and  k  blocks  is  the  nusibar  of  ways  that  n 
aambere  can  be  aurranged  (nl)  raised  to  the  power  k,  the  uioAier  of 
blocks.  These  (n!)  values  of  W  contsdn  aultlpllcltles  because 
many  of  the  shufflings  will  give  rise  to  what  Is  equivalent  to  Inter¬ 
changing  the  colusms  of  numbers  In  Table  1.1(a).  The  number  of  such 
different  Interchanges  is  nJ.  Thus  the  nunOier  of  distinct  equally 
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likely  value*  of  W  1*  et  noet  (nJ}Va'.  or  (n!)*^  At  moBt  refers 
to  acoidentaJ.  equalities  of  W  duo  to  the  date  being  rounded.  An 

exact  test  for  the  5'iH  hypotSsseis  usuid  be  t-o  ccs^are  tile  nctual  obscr” 
ved  YoJue  of  W  frou  the  eiperlnent  with  aill  the  other  values  of  sf 
froo  the  ih’iffllog*.  If  not  more  thaui  soae  pre-assigncd 

of  these  (n!)  values  of  W  exceed  the  experlaenteJ.  value  ooe  re¬ 
jects  the  null  taypotbeele  at  the  aft  level.  To  avoid  the  computation 
reguliad  to  find  all  these  valuee  of  W,  Pttssui  Inveeti gated  the  dis¬ 
tribution  of  W  vhea  the  underlying  distribution  deviates  froe  normllty. 
These  are  two  different  d'' stribsitions  of  V  being  dlseusoed.  One  is 
the  distribution  of  W  ubicb  Is  associated  with  tho  underlying  diatrl- 
butlon  and  hence  dees  not  depend  on  the  Individual  values  of  the  sa«q>le. 
I;h9  ether  is  the  distribution  of  W  obtained  frees  the  permutations  and 
depends  directly  on  tho  samsile.  ISjase  dletrlbutiooe  will  be  denoted 
by  unconditicmal  and  ccmditlonal  respectively. 

Under  the  assuoption  that  the  underlying  dlstributloo  is  noraal, 
the  uaoondltlcsival  ddstributioti  of  V  is  a  beta  with  (n-l)  and 
(k-l)(o-l)  degrees  of  freedoa.  If  W  dsvlated  very  little  from  the 
beta  distrlbutliw  for  any  reasonable  taderlylng  diatributlon^  then  ooe 
i..;uld  *assuBS  notmLUty”  and  proceed  to  uss  the  T  test  without  too 
iHich  concern,  slnoe  a  tc»t*  would.  In  general,  reject  reelly 

true  null  hypotheses  approxlwtcly  oi  of  the  tiae. 

It  f»a  shown  by  Eden  and  Yatew  [3]  that  there  was  good  agrsaaeDt 


between  tbs  beta  distrlbutisw!  and  the  conditional  distribution  of  W 
in  a  particular  saa^pllng  asperiaeDt.  Pltaan  [12]  obtained  the  first 


four  periEUtatlon  oonnents  of  W  ia  t.'jr!nB  of  the  flrat  four  EiomeDts 
*71,1  of  the  choer*.~tlos£ .  W’elCu  Is  tbs  i»fmr  that 


pansllel*  Pltaao'*,  obtained  only  the  first  two  ocaaente.  By  looking 
at  the  aoaentaj  PltDOn  vas  able  to  draw  certain  conclusions  about  the 
distribution  of  W  as  tl^e  underlying  distribution  deviates  froo 
norsiaXlty.  Pltuan's  aciBeDto  cen  be  Dade  to  look  simpler  and  hence 
easier  to  study  by  Introducing  a  change  of  notation.  Let 


2  r  /■ 


vhere 


3J 


-  1 


^diere 


,  i£i<*ii  - 

“kj  “  '  k~- 1 


l=l'*U.^  -  3a; 


The  odtatlOD  £"  la  used  to  denote  eunnatlon  over  all  subscripts 
contained  inalde  the  eunatloo  ooltting  all  terms  of  the  suoaatlon  for 
vblcb  the  eubecripts  are  s^ual.  Further  ’’avep"  Is  used  to  indicate 
that  the  average  Is  being  taken  over  penautations  v.-ithin  blocks.  ulth 
these  notatlonal  changes  Pltuan's  atoments  nov  appear  as 

avep(w)  -  i 


_  t  / f . .t  \2v 


k 

r* 


i?(n-l)  ‘‘ 


/r  2.2 


3  r  0^0 V,  k(n-i )(«-£)  r 

■  y 


-3, 


IV 


-  a 

X  P 


k— !Efa!K.!a 


aVL-p({W  -  avq>(w))3] 


-  8 


a'V‘ep((V  )*^) 


(£’■  oV')^' 

.  A  P_.3 _ 


W) 


& 

r." 


2  2 
c  a 

D  <J 


1*8 


+ 


rWr"s  ^  6(n-l)(n-2)(D-.3)  f'SpC,°X 


k 


96(n-2)  r  3; 


Tbeae  fonnuLit  cun  be  aisjjlll'ied  Bcanewhat  by  letting 


Also  Instead  cf  giving  the  first  four  laonents  of  W 
cumulants  of  kU.  ibese  are  tile  average,  variance, 
tloo  (or  Irurtosls). 


we  will  write  down 
skevnenn  and  elonga- 


avep(kU) 

v&rpikM) 

•ksp(kW] 

elop(kW) 


1  *  ?  P 

~  E”  T  T 

n-1  i  p  (i 


«vep((kW  -  avep(KW}  )  ) 

avep((kW  -aveplkWj)3)  .  _«  |vW  > 


evep((kW  -  avepCkW) )  } 


(a-ir 

3(avap(b:W  -  avep(kw)  )^)  )^ 


-lift 


T  T  T  T  + 


(n-i)^(n+l)  ^  'P  **  in-1)^  ^  V  <1  r  a 

^  0  0  T V  .  r«  S  t^t3  I 

n5(n+l)  1  P  <l  P  q  1  P  q  1  » 


Ir 

E' 

1 
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Of  these  four  cusiularits  of  W  --  average,  variance,  sliewnosB  and  elon¬ 
gation  --  only  the  average  la  Independent  of  the  particular  observatlona 
from  the  eiq>eriment. 

The  nean  and  variance  of  a  beta  distribution  with  (n-1)  and 

(k-l)(n-l)  degrees  of  freedoo  ar-*  i  and  -  reepectl"'ely . 

*  k"(kn  -  k+?) 

The  average  value  of  W  obtained  froia  rsundooiitatlon  Is  thus  the  eaaie 
as  the  average  value  for  the  beta  distribution  with  the  proper  number 
of  degrees  of  freedom.  When  the  variance  of  the  beta  distribution  If 
e({uated  to  the  vsirlance  of  W  one  gets  the  condition 


2  Cm  2  2 
■  \  E  T  T 

n-1  1  p  q 


2(k-l) 
kc  -  kt2 


2  2  1 

When  all  tu?  block  variances  are  equaJ.  then  t  »  t  •  -  for  1  <  p, 

P  4  »  - 

q  <  k.  Itius  the  left  hand  side  for  equal  block  variances  beccaaes 

2  _  2(k-l) 

(n-l)  k^  k(D.l)’ 

U>e  right  hand  side  of  the  equation  can  be  written  as  g 

which  Is  always  ssaller  than  but  by  an  amount  which  becomes 


sbmXI  as  k  and  n  Increase.  The  value  of  the  left  hand  side  de- 
2(k-l) 

creases  from  to  ler  j  as  one  block  variance  bscosaes  much  larger 

than  the  others.  Thuc  the  variance  ot  V  Is  not  greater  than 


k3(D-i) 


and  It  t&kea  this  y»lue  vben  the  block  vorloncos  arc  cqu^.  when  the 
variance  of  W  is  near  a  maximum,  the  coniiitlonal  variance  agrees  well 
with  the  variance  cf  the  appropriate  i^eta  distribution. 

If  a  few  of  the  block  variances  sum  very  large  relative  to  the 
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others,  then  the  conditional  variance  of  W  wtU.  be  smaller  than  the 
variance  of  the  corresponding  beta  distribution.  Pitman  sugg^ste  three 
possibilities  vhen  this  arises : 

(1^  Discard  the  blocks  which  have  large  variance  relative  to  the  other 
blocks . 

(il)  Fit  a  beta  distribution  by  the  use  of  the  first  two  momente  of  W 
(investigated  b"  Welch  (1^]).  This  will  mean  a  decrease  In  the 
Dusd>er  of  degress  of  freedoa  in  the  bet*  distribution  which  applies 
vmder  ncmllty, 

(Hi)  Make  all  block  varlaneea  equal.  This  of  course  requires  the  cal¬ 
culation  of  each  of  the  block  variances  and  the  adjustaient  of  the 
observations  In  each  block  to  nsike  the  necessary  changes  In  the 
block  variances. 

One  should  not  consider  doing  (i),  (ll),  or  (dl)  if  suitable 
noraallty  aaeuagptlons  can  be  aade,  even  If  the  block  variances  seea  to 
be  very  ’inequal. 

Sbevhart  c~<l  Winters  (9]  tested  Student's  t  distribution  for 
the  two-sa^le  problem  by  tsdilng  saaples  frcss  nomal,  rectsuigular  and 
right -trlangudar  distributions  and  found  that  there  was  only  satisfac¬ 
tory  agreesmnt  for  the  case  of  the  sample  from  the  norsml  population. 
Pearson  [7]  did  a  similar  andysts  for  sasq^Ies  drawn  from  several  nan- 
soraal  Peirscssias  dlstrlbutluii*  but  his  resuita  were  somewhat  Inconclu¬ 
sive.  Rider  [8]  has  treated  the  case  of  sejnples  from  U-ehaped  dis¬ 
tributions  and  Baker  (1)  has  siaillariy  treated  theoreticel  non -homogeneous 
populutlons . 
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2 

If  the  block  varlMceB  are  all  equal,  or  approxitately  eqjjal, 

the  second  mcaent  of  W  under  rendoml  nation  theory  will  be  too  large, 
but  this  has  an  appreciable  effect  only  when  k  and  n  are  aiaakll. 

Ti'e  vairlaQce  of  W  la  fairly  lnaeDBitl'''e  to  cbangea  Id  the  vaXuea  of 
the  block  yariaccea  when  k  la  large.  Pltoan  recosmeoda  that  the 
eecood  BaaeDt  of  W  be  calculated  and  coiqpared  with  the  variance  of 
the  aqpproprlate  beta  dlatributlon  when  k  or  n  la  leas  than 

Pltnao  fiirther  reaarka  that  If  a  beta  dlatributlon  It  fitted  by  aesne 
of  the  first  two  aKsaeDts  of  W,  the  third  aod  fourth  aoaeDta  a.**  likely 
to  agree  ueU.  provided  that  the  aecood  aonent  la  ot  too  aaall.  for 
equal  block  variances,  fitting  a  beta  dlatrlbutloo  in  this  aauiner  gives 
a  good  approxloatioiQ  prorldad  k(o>l)  la  not  too  sBeill. 

3 .  Sowe  yurther  Bemexka 

Besides  these  results  of  Pltaan,  other  vrltera  have  lovestigated 
the  anadyale  of  variance  procedure  when  the  underlying  distribution  Is 
non  -Domal . 

It  is  often  suggested  that  before  making  an  analysis  of  variance 
test  It  Is  vise  first  to  oake  a  test  for  hoisogenelty  of  variances  (e.g. 
Bartlett's  test).  Box  [8]  shows  that  when  little  Is  known  of  the  form 
of  the  underlying  distributioo,  madting  such  a  test  is  often  more  dan¬ 
gerous  than  oailctlng  it.  box  further  states  that  whan  there  'ireqiisl 
treatsaent  sizes  k^  and  the  variances  from  treatment  to  treatment  are 
suspected  of  being  unequal,  then  it  seeos  that  the  usual  arisiysls  of 
variance  procedure  that  says  to  pool  the  wichin  treatment  eetiiates 


shouli.’  be  replaced  by  a  criterion  that  tcjiloys  a  veighted  Sitlnate  of 
bt  jaen  treatment  vaxlaoce  aa  euggeated  ty  Welch  [IJJ  ajd  by  James  [6]. 

?o  lUuatrate  this,  aeaume  a  altaatl&u  In  uhlsh  vs  have  n  treatments j 
-  h  2 

the  ‘  treatment  having  onlta.  If  la  the  setlmate  of 

wviaace  for  the  1***  treatment,  la  the  mean  of  the  treat- 

n  k  kj 

meat  and  overall  mean,  then  the 

&  2  2 
criterion  would  be  ^  where  Wj^  »  ,  lo  place  of 

•5  2  2 

th*  u»i4)U.  pooled  ©stisa*.*^  criterion  “  x  )  /s  vbere 

^  n  2  ° 

•  -  -  1)*^  .  In  fact  Box  auggests  that  the  criterion 

baaed  on  a  weighted  eettmate  of  the  between  treatment  variance  night 

even  be  used  for  the  case  of  equal  grovgjs  where  heterogeneity  of  variance 

night  occur.  The  criterion  would  be  -  x  )V’j  Instead  of 

n  2  2 

the  uBual  ~  ^  because  it  see.as  isore  reasonable  to  weigh 

those  treatsmnt  nsane  which  have  large  within  treatnent  varlauices  less 
than  those  treatments  which  have  small  within  vreatment  vsLrlancee.  ’nils 
criterion  asema  sensible  to  use  if  one  first  computed  the  covariance 
between  (x  -  x  )''  and  e  and  found  it  to  be  positive,  as,  it  Is 
hoped,  would  usually  be  the  case. 
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Ch*pter  2 

DXADB  JUIE  DttDICS  FOK  TDC  MULTlDIWOiaiOIlAL  PROBI^EH 

2.1  Tvo  Treatwpt  Multiaiaen«loii«LX  Ca»e 

Ccanlilar  •  tve^my  l«yout  m  in  C1na$ter  1  vlth  k  Blocka  but 
with  QtUjr  tvo  treatMDt*  p«r  block.  Suppoce  nov  that  Mwml  cbarac- 
tarlsttcn  an  msoaund  cn  eacb  Individual,  Inataad  of  but  one,  m  van 
tba  oaM  vltb  tha  vork  of  Fltaan  and  Wslch.  tba  two'vay  claaalflcatloo 
can  ba  diaplayad  '.a  follova : 

TRBA1M9ITS 
BLOOKa  1 

i 


"•12 


22 


i  V  \2 

waan  1  **  !>  *  •  *,  kj  J  «  1,  2  le  the  ■eaau'aaeot  of  tbs 

reapooaa  of  the  treataent  applied  to  an  eleaeat  In  the  block 

and  Id,  In  general,  a  m-diaenaional  vector.  The  coiponent  of 

Xjj  la  the  response  of  the  k^  character  neaaurel  on  the  1'*'*'  block 
subject  to  the  treateent.  Detailed  con  si  deration  here  la  llinlted 

to  cel;  tvo  treatjurots ,  but  the  oethods  used  should  generalize  tc  inol^sde 
£X>re  than  two  treataients.  Thus  It  (ulgbt  be  posslbls  to  find  reoultt 
for  MUlttdlisenalonal  analysis  of  variance  that  are  suitably  analogous 
to  Pitman’s  results  In  the  univariate  case. 


.  lU  - 


2 .2  MvLltldlaenalon&l  Matheaatlcttl  Tools 

One  problem  that  Inmdletely  arlsea  in  the  multidimensional  caoe 

is  that  of  underatandlng  the  meaning  of  what,  In  a  single  dlmenaion, 

was  a  sum  of  squares.  For  exeuiqile,  ve  have  to  be  able  to  generalize 

k  n  2 

the  meaning  of  such  things  as  *1J  *  *  when  the  and 

X  are  m-dleMnslcnal  vectors.  It  Is  also  necessary  to  define  them 
In  such  a  way  that  they  sure  converlent  to  use.  One  wny  of  defining 
stau  of  squares  of  vectors  was  treated  by  Tukey  [13].  Re  makea  use 
of  qusmtltles  called  (by  Olbbs  (k])  dyads  and  dyadlcs.  Dyads  are 
formed  by  the  "product"  of  vectors  In  a  way  best  shown  by  the  following 
example : 


(a^,  b^).(a^,  bg) 


Ilhe  SUB  of  two  or  more  dyads  Is  called  a  dyadic.  This  addition  to 
form  a  dyadic  Is  performed  coaq>onentvl st; : 


b  c 

1  ^  i 

*  ■‘2 

'’l  ^  ”2  '■•1  ''b-t 

r 

e  f  ' 

'l  n 

t  k  *2  ‘2  ‘ 

F 

*  ^ 

«1  -  =2  ^  ^2 j 

«i 

ij 

k  ^2  iJ 

(e 

O 

t  Bg 

hj^  +  hg  T  Ug/ 

nieae  operations  need  not  of  course  be  restricted  to  two  or  three 
dimensional  vectors.  For  an  analyst s  of  varisons  oitK  vectors  instced 
of  scalars,  we  may  replace  the  sum  of  squares  by  s  dyadic  of  sume  of 

squares  and  cross-products .  The  suma  of  squares  resulting  from  a  set 

2 

of  k  m-dlsiensionai  vectors  will  be  a  'dyadic  with  a  conponents. 
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Ab  an  exaagjle  conalder  a  dealgn  consisting  of  only  two  experlnent# 
with  two  Bsasureaents  being  taken  In  each  experiment.  Suppose  that  the 
results  are  as  foil  owe : 

gxpertswnt 

Character  I  H 

A  5  3 

B  3  1 

nia  d^jrad  idUeh  Is  the  contribution  to  the  total  sums  of  SviUares  and 
cross  products  from  the  first  experlsient  Is  the  product  of  the  vector 
(5»  3)  with  Itself.  Thus  ths  total  sums-of-siiuares'and-erosB-products 
8  Is  the  dyadic  with  four  eoigsonents 


S 


(5,3).(5,3)  ♦  (2,l).(2,l) 


r  ”Vf‘  ')-r 

\15  9/  \2  1/  U7  10 


2.3  Dyads  aa  Applied  to  a  2  x  k  ftcperlaent  In  ■  Dimensions 

tie  will  now  proceed  to  apply  these  rules  of  operatloo  of  dyads 
and  dyadlea  to  the  2k  vectors  given  In  section  2.1.  let  c^  be 
the  difference  vector  F>»rtb»rs»re  let  5^,=^  "or 

1  <  1  <  k  be  the  set  of  dlffereoce  vectors  when  the  x, j  are  permuted 
within  block  1.  In  this  way,  «  +1  If  the  difference  vector  is  as 

found  frfim  the  O  JServed  data  and  -1  If  the  dllference  vector  is  se 
found  after  Interchanging  the  Tf-ctcrc  in  block  1.  Tuo  t-.  '  ai  "sua  of 
squares"  of  the  dlffere'iceB  for  such  au)  arrangement  can  be  partitioned 


as  follows : 


Let 


A 


k 


"  HKSmiAL 

k  1«1  1  1 

I  *  ? 

,  aisa  C  -  RESrrUAL 


Vurther  let  us  denote  the  con^jcmente  of  the  n-dlinensloaal  vector 

2 

*1*12'  ’•*•  *i*i«^‘  A  Is  e  dyadic  with  n  term* 

obtained  by  adding  tc  dyads,  the  1^^  one  of  vhlcb  is 


Ttuis 


A 


^il 

‘il*12 

■  •  •  “ll“la  ’ 

2 

‘ii^ia 

*12 

*12“la 

; 

‘ 

* 

• 

•  / 

‘ll*l* 

•l2*iB 

•L  / 

1-1 11 


1^*11*12 


lSl*il*12 


E  2 

lil‘l2 


”  l-l^ll^la’ 


iSl®!?*!” 


k 

^  OL  a 
i.sl  11  iB* 


k 

l^*12*ls 


r  2 
lil*iB 


/ 


B,  oti  taa  otter  hand.  Is  a  dyad  which  Is  ’foraed  by  taking  ^  tlases 
the  product  of  a  vector  D  by  itself  where  D  le  obtained  by  Bujoalng 


the  »-dlii»nslo.-s8.1  vectors  e^c. 


over  the  k  blocks . 
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k  k  k 


<i5iVii>{i5i‘i'i2> 


B  -  ^  D.D 


2 


k  k 


l\m\ 


/ 


Pltmn  conalderadj  In  tha  univariate  case,  wbat  vould,  in  the 
present  notation,  be  denoted  by  b/(B  +  C)  where  each  would 

then  be  a  scalar.  Since  B  -s  C  A,  the  ratio  nay  also  be  written 

B/A.  This  ratio  In  the  uolveu-late  case  Is  merely  a  (dlsaDSlonless) 
numerical  quantity.  The  "ratio"  which  will  be  used  here  in  the  raulti- 
dlasoslonsl  situation  is  B/A  looked  at  In  such  a  way  tliat  A  remains 
constsmt.  A  convenient  constant  for  A  is  the  m  x  d  dyadic  whose 
dlaabnal  elements  are  all  one  and  whose  off  diag-onsl  eleaents  arc  all 
zero.  A  further  discussion  of  this  “ratio"  B/A  Is  given  in  Chapter 

3.  It  would  then  be  enough  to  compute  moments  of  this  new  dyad. 

Purtiiennore  one  would  not  need  to  attimpt  to  attach  a  aeaning  to  the 
ratio'  or  an  arbitsnry  dyad  and  an  arbitrary  dyadic. 

llils  aiB^illflcatlon  requires  that  r.h,  ,,  1  »  1,  2,  k  and 

k  „ 

J  «  1,  2,  ....  0  ce  ren^aced  by  new  numbers  t,a'  such  that  «  1 

j,  1  'J  Ip 

for  1  <  p  <  ui  and  -  0  for  1  <  p  /  q  <  o.  Toe  class  of 


^or 


transf jraBtioiiB  whicii  aseet  thees  requlrerafinte  is  a  linear  one. 
the  special  case  of  two  Ali&ension#  a  transformation  satlsC.v i .;g  these 

requirements  Is  (jlven  by 


and 


Li/<‘ 

J'  \ 


1; 

I,  a  a 

I  2 

P-1  pi 


11 


q2 


£.1 


inila  la  only  one  of  many  transformations  such  that 


3 

t  *'* 

P-1  pi 


£  ..2 


p»i%2 


1  and  2  a'  a'_ 
pel  pi  p2 


0.  Ifce  class  of  transforma- 
tioaa,  of  which  the  transfonaatlrai  given  by  equations  2.1  Is  a  member, 


eamalsta  of  all  a 


■pi’  “p2 

a*.  ■  I 

pi 


auch  that 


COB  0  -  a'  sin  0 
pi  p2 


2.2 


p2 


a'  sin  9  *  a'-  cos  9 
P2  p2 


vtaare  0  <0  <2%.  Mathesiatlcal  definitions  for  the  claaa  of  trans¬ 
formations  for  the  ^eoersl  a-dlrnsnelonal  situation  follow  directly 
froo  the  two  dimensional  definition  given  by  2.1  and  2.2. 


Chapter  3 

uiLicur»3ioff  or  Tnis  criti£kioh  b/a 

3.1  rt .1 1 1  or. irig  of  the  of 

Consider  for  the  present  the  general  univariate  case  of  a 
treatroente  in  k  blocks  as  dlecujseed  by  Pitman.  We  ai*e  given 
the  measurements 


*12 

*ln 

*21 

*22 

*2n 

■Ta 

*k2 

*kn 

vbere  x.  .  la  the  observed  response  of  tho*  blork  subjected  to 

i  J 

the  treatment.  'Hie  most  general  tjuedratlo  fora  relating  these 

Xj^j'e  la  given  by 

h  n  k  n 

«  “  iSl  j5l  1^1  A  ^IJIJ^IJ^IJ  3.1 

In  order  for  q  to  be  invariant  under  Interchanges  (shuffling)  of 
Blocks  among  themselves  and  of  treatments  among  tbemseives,  certain 
ccaidltlons  involving  the  a,  are  necessary.  To  detonalne  these 
conditions,  suppose  that  JJ  is  fixed.  All  tersM  in  q  for  which 
this  is  true  are 

k  k  k  k 

i5i  i5i“ijij*ij’‘ij  “  *  i,?-i“iji/i/ij  3.2 

Suopoae  all  but  two  of  :he  are  zero.  Further,  suppose  that  tiia 

two  non -zero  x^.'s  ajpear  In  the  aome  block.  Let  us  denote  these 


^*0 


y.  by  X  ,  uni  x  The  flrat  suiiaaitlon  I’et  it  be  heriuted  by 

IJ  pj  pJ 

£,  )  on  the  rl<^it  hand  side  of  eoneffor';  'i-'~  Hi’i  have  cy.a..'tiy  or.. 


non -zero  term  a  ,  ,x  .x  ..  tJuiipose  nov  the  blochs  are  i.lnU’J'led  no 
pjpj  pJ  po 

that  block  p  no  longer  contains  the  two  non-cero  'i'he  value 

i  J 

of  tlus  t»ily  non-zero  term  In  E,  will  now  be  a  ,  ,  , -X  ,  .x  ,  ,  where 

1  p  Jr  J  p  J  p  J 

block  p'  la  now  the  block  containing  the  non-zero  Unleas 

C  will  not  remain  unclianged  under  shuffling 

pJpJ  p'Jn'o  JJ' 

of  the  block! .  If  we  now  euppoee  again  tliat  only  two  of  the  ^ j  j  '  * 
are  non-tero  but  this  time  assume  that  those  two  non-zero  x  'a  are 

J-O 

In  two  separate  blocks,  the  second  summation  (denoted  by  Z^)  on  the 
rlijht  hand  aide  of  equation  'i,2  will  have  exactly  two  con-zaro  tem.s 

shuffled.  It  is  necessaj-y  that  the  coefficients  a  ,  ,  and  a  ,  ,  be 
'  PJqJ  kJpJ 

constant  for  fixed  JJ.  ISius  it  Is  necessary  that  equation  3.2  bo 
of  the  form 


k  k  k  k  _ 

iSl  r5l'^lJIJ*lJ*lJ  ■  ®JJ  iSl'^ijij  *  1,F.1^1J*IJ  3- 

By  reversing  the  roles  of  the  blocks  end  treatments  and  roj -eating  the 

above  argument  it  follows  that  a  general  quadratic  form  satisfying 

the  required  sjnmaetry  with  respect  to  blocks  and  treatments  can  be 

partitioned  as  follows: 


k  ii  .J  n  k 

i  f,  t  b.f,  S'',x. 

1=1  J-,i  IJ  - 


'IJ 


lii  J 


,5>i’'i 


j^i 


k  r. 

S,i=i 


j'-ij 


3 --i  The  Adeq'jacy  of  Studying  P/J  ,  provided _ A _ 1.  s  scaled  as  I  ■ 

From  the  preeedlog  argiwients  on  the  part.ltlorii,ng  of  »  quailratic 
form  it  follovu  that  the  total  eua  ol  8qiJare:i  ol  the  x,  '  a  can  be 
partitioned  into  four  parte.  With  some  algebraic  manipulation  the 
partitioning  can  bo  reduced  to  a  suai  of  squares  involving  the  mean  of 
the  X  'aj  one  involving  the  nteans  of  the  treatiaente  (x  ),  one 

1  J  •  J 

Involving  the  means  of  the  blocks  (x^  )  and  finally  one  Involving 
both  the  treataent  means  and  the  block  means  (residual  or  error  a'ia 
of  squares  J.  Without  any  prior  Infonnatlon  on  the  mean  of  the  under- 
ij-ing  population,  any  criterion  which  tests  for  the  difference  between 
treatoento  should  be  independent  of  the  overall  loesui.  Since,  moreover, 
rsuidcanlzattoos  within  blocks  leave  the  block  sun  of  squares  imchanged, 
we  need  concern  ourselves  only  with  a  r;rlterlon  involving  the  treatisent 
sum  of  squares  and  the  error  sum  of  aqioare.s.  These  two  subib  of  squares 
were  denoted  by  B  aid  C  in  Chapter  2.  A  ratio  Involving  B  and 
C  is  chosen  arbitrarily,  following  the  Pltmarj  approach.  The  most 
genereJ.  ratio  Involving  P.  and  C  la  of  the  i'orm  (cffl  +  PG)/(7B  +  6C ) . 

Consider  for  the  puipose  of  i  j  >  .istrsti  on  the  “peclal  case  where 
only  two  treatments  ai-e  under  test  and  where  only  a  single  response  is 


being  in>>ttsd;-ed.  buppoue  that  we  scale  the  oeabureaects  such  that 
B  +  C  1.  We  vtil  lierote  H  by  and  C  by  Cj  uniler  this 

\  na  - 


Sv*x>pCBe  an  fixperiiiscrit.  .joris of  blocks  with  two  treatnicntti 

per  L)locx  ajiri  one  ...perlicent  on  each  t re&t.'.-ierit  In  each  block.  7*A»e 
reiipoiise&  can  bo  exoreo^ed  by  five  nuaiberB  which  are  the  difrerences 


of  the  ytexda  In  the  five  blocks. 


Kurthencore ,  these  five  differences 


call  then  ,  cen  be  displayed  as  points  on  a  line  as  shovn  in  Figure 
1.1.  If  we  nov  change  all  the  signs  of  the  a^,  ve  obtain  five  more 

- sp-i - 

Figure  3-1 

points  to  place  on  the  line  which  are  the  lame  dietsmce  fron  0  In 
the  opposite  direction  aa  shown  by  Figure  3-^>  By  a  rescaling  of 


t;; — J - fpf; 

Figure  3.2 


these  differences  a^  ve  can  find  a  unit  of  oeasui-esient  such  that 


5 

E  af 

1«1  1 


1. 


When  this  has  been  done,  any  paurtlcular  randooizatlon  R  giving  rise 

to  differences  *^*2>  *3*"^' 

or  -1  (depending  on  whether  occurs  without  with  a  change  of 
sign)  Is  such  that  squares 

la  B  *  C  and  any  criterion  based  on  the  a^  other  than  one  which 
has  B  +  C  “  i  (or  sosie  other  convenient  constant)  or  a  m'Htlple  of 
B  +  C  In  the  denoHilnator  will  have  a  different  denossinator  for  each 

randosiJ  tttlon  chosen.  The  numerator  can  have  any  linear  combination 
X  ’ 

and  “  i?" 


of  1 


A  reasonable  criterion  Is  thus  of  the 
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form  (aB  +  6C)/(B  +  C).  However,  by  finding  the  randmni tation 
cunulants  for  «  special  case  of  this  criterion,  nauaely  3/  (B  +  C)  or 
B/A  scaled  such  that  A  •  1,  one  can  then  find  the  randainlutlon 
cwnilsuits  for  the  itore  general  criterion  (oB  +  pc)/k.  The  formulas 
relating  the  rauidOBlzatlon  cunulante  for  the  general  critarlss  (aB  -f  ^C)/A 
in  terms  of  the  cumulacts  for  b/A  are  now  derived. 
avep((3B  +  0C)  «  a  avep(B)  t  p  avep{C) 

•  a  avep(B)  +  p  avep{3  -  B) 

-  (a  -  P)  avep{B)  +  p 
varp(aB  +  pc)  ■  varp((a  -  p)B) 

■  (a  -  p)®  varp(B) 

8kep(a®  +  pc)  -  Bkep({a  -  p)B) 

-  (a  -  p)^  ekeplB) 
elop(aB  pc)  m  elop{(a  -  P)b) 

-  (a  -  p)**  elop(B) 


3*3  Multlvarlete  Analog  of  B/A 

In  the  multrresponse  situation  the  condition  B  t  C  ■  1  becosses 


B  +  C  «  I  where,  for  the  caee  of  a  responses,  I  Is  a  m  x  a  dyad 


with  I's  on  the  diagonal  and  O'a  elsewhere.  Again  it  will  make 
DO  difference  whether  we  are  cocalderlng  Bj  (b  under  the  conditloo 
B  +  C  »  l)  or  Cj.  or  sooe  linear  combination  of  Bj  and  C.,..  As 
waa  the  case  i.n  the  univariate  situation,  the  randosalzatlon  cumulants 


for  the  criterion  B  with  suitable  scaling  so  that  A  -  I  are  suffi¬ 


cient  for  obtaining  the  cuniul».its  of  smy  arbitrary  linear  ooml’lnation 


of  B  and  C . 


avep{Q®  +  pCj 


WflaB  *  pc) 


*ksp[CEB  +  pc) 


nsally 

alojifiXB  +  pc)  • 


>  avepjo®  +  p(i  -  B)} 

“  M*  ’»•  •»  p;rj 

-  PI  +  (a  -  p)  av*p(B} 

-  wptOB  +  p(r  -  B)3 

-  Y»rp((a  -  p)B  *  PI) 

-  vajT>((a  -  p)b) 

-  av*p{(a  -  p)B,(a  -  p)B)  ~  •v»j>((a  -  p)B)^v»p((a  -  p)B) 

-  *vep((a  .  p)2b.b)  -  (a  -  p)*  aveplBl  ^v«p{B) 

■  (a  -  p)  [av6p(B,B)  -  avcp(B}^vep(B)] 

-  (a  -  p)®  VBrp(B) 

■  ekepIOB  +  p(l  -  B)) 

«  elwpC  (a  -  p)B  +  PI) 

-  »kep[(a-p)B) 

'  avapifCQ  -  B)B  -  .v«p((a  -  p)B)J,{(a  -  p)B  -  avep{(a  -  p)3)). 

((a  .  P)b  -  avep{(a  -  p)b)  )) 

.  *«p((a  -  P)[B  -  avep(B)j.(«  -  p)[B  -  a«,p(B)  ). 

(a  -  P)[b  -  avepts)]) 

*  (.^  -  P)^avepf{B  -  avepCB]  ).(B  -  aYap(D),(B  -  avapfB)  )] 

(a  ,  p)3  stepCB) 


aIop{aB  +  P(l  -  b)) 
eiop{{a  -  p)b  +  PI) 
eiop{(a  .  p)b) 

avaF([(a^)B  -  av*p((a^>)El  ].t(Qr^J)B  -  avep(  {a.p)B)  ]. 
{(a.p)B  *  av«p((a-<s)Bjj.i(aHp)B  -  avep((a-p)Bj  ]) 

-  3  rarp((Q-P)B).YRrp((aHB)a} 
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■  avep((c!-ft)[B  -  av«p(B}  }.(a-il)[B  -  ».vep{Bl  l.(a-<»)tB  -  avap(Bl  ). 

1)  -  »»rp(B) 

■  (a-0)^  »Yop{(E  -  avcplEj  )«(B  -  iiY«p{B)  ).(B  -  evep{B)  )« 

{B  -  av«p(B)))  -  3(a-<J)*  YMp(B).vmiT>{B) 

«  (a  -  g)**  •J.op(B)  . 

tii»  derivation  ot  th«  flrtt  four  etajulants  of  B  for  all  pennuta- 
tlODS  of  Plgn  li  given  in  Chapter  H.  In  keeping  with  the  notation  of 

Chapter  2,  the  dyad  B  ie  foneed  t>y  taking  the  product  of  the  vector 

k 

of  SMM  Of  differences,  with  itself,  divided  hy  k,  the  suaber 

of  blocks.  This  can  be  vritten  as 

®  “  eC  i5i®i  )  ( i5i=i ) 

The  results  of  Chapter  4  are : 

1  k  1 

avepiB]  «  E,  c  ,c  “  r  * 

k  v«l  p  p  k 


varp(B) 

skep(B) 

el  op{B} 


-er  C*  (c  aC  aC  ,0  +  C  aC  aC  aC  ) 

i^2  p^^.i  p  q  p  q  p  q  q  p' 


1  *■ 

-5-  E"  (c  »«  *C  aO  ),(c  «C  -VC  aC  ),(c  .0  +C  aC  ) 

k3  p,,.r.l  P  <1  1  P'  •  -i  r  ^  p  p  r^ 


V'p>‘^V"r^  V 

=8*‘=r)«('e*V 

(c  aC  *  )}  -  £*'i(c  mC  ¥  C  mC  )*(c  gC  +  C  gC  )c 

tjr  r  <5.  P>  'p  ^ 

(c  ,c  +  o^.c  ),(c:  ,C  +  C  aC  )), 


invaria.;;Uij  aaicl  CuPiuilants  over  Rotation  of  B 


in  order  to  provMe  r*  proceoiur«  for  testing  whether,  in  any  given 
inetance,  the  two  treatxnenta  are  really  alike.  It  is  desirable  to  be  able 
to  compute  a  suitable  quantity  from  the  data  and  compare  it  with  son* 
predetermined  value.  '1‘he  quantity  that  Is  computed  in  each  Instauce 
shoiU.d  be  independent  of  the  particular  choice  of  coordinates  In  tertca 
of  which  the  rebpons-’-e  are  expreated.  The  e' ements  of  the  dyad  B, 
even  vita  the  reetrlotion  that  A  be  equal  to  a  given  constent,  will 
depend  on  the  choice  of  the  tranafc-rmtlon  ueed  to  aiaJ-:e  A  the  parti¬ 
cular  conetant  va  wish  (most  conveniently  3).  However  the  trace  of  B 
1b  Independent  of  thi‘  choice  of  transfr>r«vit.icn  which  maKes  A  equaJ 
to  acazte  convenient  conetant,  'i’hi©  {(.  now  nhown. 

For  convenience  l.i  notation  let  the  obeervatlone  after  scaling  to 
aalce  A  •  3  he  giver,  by  (a^,  b^)  In  block  1.  Then 


k  k 

i  iol  1' 


*  "*  ‘  / 


and 


traces  -  ^  (  jii-*,  T  +  (  )' ) 


Let  MJi  r.ov  writi-  5  lu  n  CvCiiu-inauc  KysLcm  V®  '  i  b*)  whicn 

has.  beun  rotated  ihi  oii^h  an  angle  0  froin  (a,  b).  Then 
a  »  a’  ecu  6  -  b’  aln  6 


a*  sin  +  b'  C03  d 


Hence 


trace  5  -  ®  S  +  ^|  cos  9)]' 

k  k  k  k 


k  o  k 

(i5i‘P  ^  (iVP  • 


?nbus  trace  B  doea  sot  depend  on  the  eiigl9  8  through  which  the  aixea 
were  rotated  and  la  therefore  Independent  of  the  orlentatloo  of  the 
retpooae  variablea. 

To  obtain  the  cumlanta  of  trace  B  over  perautatlona  we  need  ts 
define  eon  new  quantltiea .  let  1^0,  denote  a  linear  form  of  the 
elemente  of  B.  Here  o  la  a  double  subscript.  la  particular,  If 

Is  such  that  1^^  is  trace  B.  then  c«e  only  needs  to  express 
the  peraattatlon  cumulants  of  In  terns  of  the  penoutatloo  cumulantt 

of  B.  Denoting  the  elements  of  by  1^,  and  those 

of  by  Bj^2'  ®a' 

avepCl^B^J  -  aveplL^B^,  i  i  ♦  L^^B^g) 

sveptB^j^)  +  avop(B^g)  +  svop[Bgj^)  ■*  avep(Bgg) 

The  right-hand  i..iue  of  3.5  is  e^iuivalcnt  to  operstlag  on  aveptB)  . 

Thus  we  have 

aveptl^Sj^} 


4j(svep[B)  . 


3.6 


Slallarly 

varp{L^^^i 

-  I^L,Jvarp(E))^ 

3.7 

'-bra' 

_  T  T  r  \ 

•a  O 

elQp(lj^l^^ 

3.9 

Thus  th«  first  four  perautAtion  cu»uLaknt»  of  trace  B  can  te  (Atolned 
froa  3-6,  'i-7,  J-8  and  3-9  by  replacing  the  o^aratoro  Ijj, 

and  by  ao  operator  which  has  eleiaents  m  m  1  and 

I^-X^  -O. 


2 

K 


■  K  D  ‘  C  D 

,  i  O  k  *  w  > 

Tarp{trace  8)  ••  -W  ):  E"  (c„»c.<.c  «c,^  +  c  ,c  .c  ,o  )  :  (  ) 

1/  p.a-O.  PaP'i  P4<IP  Voi/ 


3.10 


1  On 


^  ^Vu  ^  v\^ 

k  p,q«l  ^  . 


3..U 


Btallarly 


s)wp(trace  B)  -  E“  (a  a  ♦  b  b  )(a  a  +  b  b  )(a  a  +  b  b  )  3-12 


&cd  finally* 


elopCtraoo  B)  -  E"  (a_a  t  b  b  )(a  a  +  b  b  )(a  a  +  b  b  )(a 

k  p,ii,r,8«.i  i-  'i  y  'i  y  y  h  ■* 


a 

r  s 


ig  a  ^ 

-t-  E"(aa  ■♦•bb) 

k**  n.a^l  P  »  PI 


M3 


The  ease  of  eoaiputatlon  of  the  pentiutetloo  cuesulante  of  trees  B 
higher  than  the  first  will  depend  on  our  ability  to  write  thee  in  a 


way  idiich  requires  the  Bunmlng  of  a  emali  nrunber  of  terns  rather  than 
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as  they  ore  given  Ic  3.i2  aad  3-13-  Since  trace  B  1» 

invariant  with  resp^^tr  to  the  orlentatlcc  of  the  (a,  b)  oxee  it 
eenoie  Aeslraule  to  coelute  the  cuoulont*  of  B  avero^l  over  rota- 
ticcs  of  the  coordinate  ayetea  to  vtaloh  the  vectoxa  deflnlog  B  are 
referred.  In  this  way,  Ute  eXoiwetc  af  these  erexeced  cuaulents  of 
B  vill  be  iUTorlaot  ic  the  aene  way  that  trace  B  is  luTorlant. 

It  tume  out  that  the  eleeents  of  the  evare^ied  cuauloats  of  B  can 
all  be  expressed  os  linear  coedbliiations  of  seyen  basic  iavarleats. 
nicse  lavoriaata  are: 


+  bj^) 

pel'  p  p' 


k  L  k  I. 

I",  (as  -f  b  b  )  ;  and  E"  ,  (a  b  -  b  al 
p^q-d.'  p  q  p  q'  '  Pjq-'J-  p  q  P  4 


Since 

avepCtroce  B}  >  ave  avepCtroce  B> 

e 

-  are  l^CevapCbJ  4. 

0 

m  I^[ave  avepCBl 
0 

ftod  It  oan  b«  siiovii  Id  &  elmlXftr  t^mhloa  tmX 

vorpCtrace  B)  »'  l^^Iave  vaip(Bj  j  , 

0 

skepttrace  3]  «  8kop{B]  , 

and  elcp{tia««  33  »  '5  Jopyj 

9 


then  the  oussAlat-its  of  tree*  E  cen  be  found  dirt  ctly  froia  the  pemuta- 
tloo  cuMulants  of  H  averaged  over  all  orlentat-i  as  of  tbs  andcrijlBg 
coordinate  systeBs,  In  this  tbs  cassiLsBts  of  trace  B  con  be 


ejfpreneed  as  Hriear  coabitiatlons  cf  the  seven  basic  invariants  Kiven 
above.  iUl  these  Invariants  can  be  coi^puted  quickly  and  easily  In 
any  si>ecif Ic  case . 


3.5  Parlvatloa  of  Cuaulants  of  Trace  B  lu  Teras  of  a  Set  of  Invariants 
Hie  foraulas  3'ld.,  3.12  and  3*^3  for  the  second,  third  and 
fourth  ctusuLsnts  of  trace  B  can  be  expressed  In  the  acre  eoaputa- 
tlooally  convenient  way  described  In  Section  3^‘^  applying  the 
rasults  of  Chapter  One  can  easily  write  these  ctuailants  In  terms 
of  the  set  of  Invariants  that  sire  obtsdned  for  the  four  cumulants  of 
B  averaged  over  rot.atlon. 

avepCtrace  B)  •  j 


varp{trace  B] 


2  ««/  V.  V 

‘  I  ^ Vq  '  W 

m  S-  ^(a®a^  +  2  a  b  r.  b  + 

l>CQd  T>0^ 


p  p  q  q  p  q 


d 

a  -f 

P 


3.15 


Tbsas  results  can  aJ.so  be  arrived  at  dlrectl;/  without  the  use  of  the 
-results  of  Chapter  5,  for 


E"(a  fa  +  b  b  )'■ 
i  p  q  p  r 


(1  -  r  -  2  K  a^^  ■(-  (1  -  E  b^) 

'  jP'  iPP  '  iP^ 


V  k  /  2 
a'  -  E  (a  +  b  ) 
1  P  P' 


'  3i.  - 


skcp(tre.ce  B) 


fl  k  . 

-■T  i"(a  a  +  6  fa  )(a  »  +  t  b  )(a  a  f  b  b  ) 

1^3  1  P  '1  p  q''  p  r  p  r''  1}  r  <5  r' 


-  ™  +  3a^fi.beb  +  3ababb^  +  bfb^^) 

^.3i  Pli'  "^p-jcirr  PP<iqr  V  i  r' 

■  ?<*  - 1  •  i  ■  I  ‘  i 

-  ^  *  I  '--i  H<>‘  •  i 


-  ^2  -  3  +  2  ^(a.®  ^  b^)^3 


1  V 


X  V 


3.X6 


RUi  ntult  also  can  be  ottalaed  directly  for 

+  3*^cfcab  +  3*b»bb^  +  b^^®) 

I'pqr  pqqrr  ppqqr  pqi"^ 

"1-3  -  I  -  3  -  6  £“*^  a  b  -31  aS)® 

iPi  iP  ipqq  ippqq  ipp 

-3  ^*^5*0  ■  ^  ■  3  ^  *5*0  +  i  -  3  ^  b^ 

xppq  ippqq  ipp  ipq  ip 

••  a.’-3i^a'*  +  2Ca^-3E  aS)^  v  3  E  a^b^  +  6  ^  a^b^  -  3  E  aSi® 
iP  iP  iPP  iPP  iPP  iPP 

-3  i  a^^  +  3  r  a^**  +  6  E  a^b**  -  3  E  -  3  E  b*^  +  2  £  b^ 
i  PP  iP?  iPP  iPP  iP  iP 

-  2  -  3  +  2a^?  +  b^)  +  2  ^a^  +  3a ♦  3a^^  +  b^) 

i  P  PP  P  P  P  P  P  p  V 

“2-3  +  b^}^  +  2  E(af  +  b®)3. 


1  P 


X  P 
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elopttrace  B)  »  '’ci’^s ^^*r*3* 

-  i|  £”(*  »  +  b  to  )'‘ 

^4  j  '  p  q  *.  qi 

“  ^  +  4  i.^t,^a,  to  «  to  +4  e?».  to  a  to  to^  +  2  a  to_a  to  a  to  a  to 

^**jL  P«r*  pqrraa  pqqrra  ppqqrre* 

+  4  a  to  a  to  to^^  <■  bS^^to^)  -  if  ^'(a  a  +  b  to 

ppqqyg  j^4  p  p  (J 

"  I"a^a^a^a^  -  16  ^a**a*^)  ♦  ^48  £"a^a^a  b  a  b  -  l6  E'’a^  a\  ) 

^4'  ipqr*  iP4  jj4'  i  pqrrss  x  ppqo. 

*  ^48  I'a^a  b  a  b  b®  -  l6  r'aS>^aS>’^)  +  ^448  £"a  b  a  b  a  to  a  b 
4'  X  p^<irrB  i  PPQ.^^  1,4'  j.  pp^qrrsB 

JC  & 

-  X6  lc"a^)^aS3^)  +  t448  J^'a  b  a  b  bS^  -  l6  T'a  b^a  b^) 

1  p  p  4  q'  J.4'  j.  pp^<Jir8  X  PPQ^' 

+  iw4a  •  i6  r  bV*^) 

j^4'  ipqre  X  pq' 

fVOB  Table  5»7  C>\apter  5>  one  then  obtains 

alop(trace  B)  -  ^^{(48+48)  +  (.108-24  -24-04-103)  |(ap+bp)^  +  (l2Of4a+48-4v3-l20) 

-2+6+15) 

(£(a  +>>  +  (20+14+4+8+14+20 )E“(a  a  +  to  to  )^  +  (-12-6-20+8-6-12) 

1  P  P  1  P  <l  P  r 

r'(a  to  -to  a  )**) 

X  '  p  q  P  ^ 

-  if(6  -  18  £(a‘'+b'')’='  +24  l(a‘^+b‘^)^  -12  l(a^+b^)^  +3(l(a‘’-+b'=)‘"i'-’ 
k'r  l'  p  p'  r  P  P  X  P  P  1  p  p' 

»4, 


^"(a  a  +b  b  )■*  -3  £"(a  b  -b  a  )’] 
I'P^  PQ*  PQl  * 


3.17 


A^in  ^.his  result  car;  be  obtained  sliructl/. 
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K  2  2  ?  2 

r'*a*"a  a  ~a 
i  P  r  ti 


,  ^  ^  6  ,  J.  fi  'ivr; 

L-o2^a  1-  o  t  Q  -  6La  +  3(£a  ) 

1  P  1  P  X  P  - '  p** 


^2?  .  62  li  22  22.  ^2  ^  \ 

Z"^  n'a.  1»  a  b  -  -b  i;  &°b  <•  I  a  E  a^  -  E  a  h  *  h  t  ^  ^  2  E"i^s  a-^D 

ipnrrss  iPP  ipipp  iPP  iPP  1PP«*1 

r"a^ft  b  a  b  b^-  -  if  aSs^  +  2  E  «S)^  ♦  (Ea^^)^  -  4  E  +  2  aa®b**  +  ^"a\  a  b'^ 

iPP  iPP  'iP  P*  1l  VP  X^P  P  X  P  P  (I  ^ 

-  -  6  E  +  3(1  A?)^ 

iPP447raa  iPP  i  p  p 

E"a  b  a  b  -  -  b  T  a^^  +  E  aS^  fb*^  -  E  a^^  +  E  aS»^  +  2  E-a^  a^ 
iPPlir*  iPP  XPPlP  iPP  XPP  iPP4m 

E*b^^^Df  .  1  .  6  E  b^  +  8  f  b^  -  6  E  bf  +  3{E  bS^ 


1  P  4  r  • 


1  P 


1  P  X  P  ''  P' 

Wb*D  tbeae  fontixlaa  are  substltutad  in  the  foraula  for  elop(tra.ce  F)] 
and  the  tonus  ajr'e  soUocted  the  result  checks  with  eijiation  3>i7> 


DiRIVATION  or  CJMtJUUrrS  FOR  B  yilDfB  RABDOMrXATIOH 


4.1  Some  reiaarks  cm  Motatlon 

The  ^uaatlty  for  which  Pitaaa  found  the  moBients  corresponds  to 
vh«t  is  celled  ^  In  Chaptex-  2  In  the  nultldtmensloneJ.  situation. 
This  ‘ratio'  reduces  to  B  Itself  when  considered  with  lespect  to 
a  suitable  class  of  tranafonsatlons  with  K  ^  1,  as  deecrlbed  In 
Chapter  2.  In  the  two  dlstenslonal  case  this  can  be  written,  with 
respect  to  a  particular  coordinate  system,  in  the  fornj 


When  ttie  treataients  are  randomlied  over  blocks  the  dyadic  B  takes 
the  form 


iri-i  1  t  'i-i  1  i’ 


Before  proceeding  with  the  cufliiuants  averaged  over  raodoaiiiatione 
one  needs  nooxi  further  knowledge  on  bow  to  work  with  dyads  amd  dyadlce 
It  1b  neceasarj,  ft;  'iaaliog  with  mcjseiits  hlght_r  than  the  first,  to 
find  everagsh  uf  'jK-.vers’  higher  then  t  le  first,  of  dyads.  It  ie 
necesbary  tv-j  define  the  product  of  two  oi'  iriore  dyads,  and  honce  th** 
ss^ufti-e  and  higher  powers  of  dyads.  iituoe  a  dyad  was  formed  by  taking 


-  - 


the  product  of  two  vectors,  the  aiiuare  of  a  dyad  or  product  of  two 
dyada  results  In  a  four-dlinenelonal  array  which  can  be  displayed  as 
fnlJ.niVS  ; 


( 


a  b  \  /  e  f  \ 

a« 

SLf 

be 

bf 

«« 

ah  ' 

bh 

ca 

Cf  \ 

\  , 

dft 

di' 

( 

\ 

eg 

ch  ^  dg 

dh 

Thase  four  S  x  2  arrays  really  should  be  thought  of  as  a  single 
2  z  2  X  2  X  2  array.  One  might  vlsh  to  think  of  the  eleoents  of 
such  an  array  as  forolng  a  k-dlmensiODal  cube  with  the  elements 
theamelves  as  vertices.  In  the  sajne  way  hi^er  powers,  or  more 
esneraily  aiultlple  products,  can  be  represented.  J^urther  notatlooa 
will  be  defined  as  the  confutation  of  cuzailants  proceeds. 


1».2  CuBBilants  of  B  utide  ■  Rapdoml nation 

llie  first  cunKiiant  of  B  for  the  m-dliaenslfwjai  sltuatloc  under 
randonlzatlon  is 


ave 

(permutations 
within  blocks)] 


£^lil*l“ll^^l-l‘^l“l2^  S^i5l*l*12^ 


,  k 


^lSl"l‘ll>ft5l«i*ia)  U 


r:{,iE,€,a  _)(  E.t.a,  ) 
k'ial  1  l£''i>«l  i 


P 


It  1.S  uijderstoofl  ttist  1 -=  the  reattilnder  of  this  chapter  "avep*  luflles 
averatiC^e  over  ons  vlt-Mfi  l5ic?ck.«. 
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To  make  tJile  expression  uiore  aanagcable,  let  us  write  It  In  vector 
notation  by  letting  c,  represent  (a.^,  a^.^,  Then 


avep{B) 


Recalling  now  hov 

avep{«^€^ 


^  k  k 

k  k 

"  k  li  j£l‘^l*“^  *''ep(c^ej) 
uao  originally  defined,  It  Is  clear  that 
)  -  1 


avep(«^€j)  »  i{+l)  +  H-i)  -0  for  1  ^  J 


k  k  k  ^  p 

Making  use  of  r"  again  to  nean 
vs  (^taln 

k 

aven  t.c.c..c.  +  ^  .  ?''.av«nff  ^  . 

J  t  J 

k 


ayep(B)  -  ^iSi*vep  c^c^c^.c^  y  j;;^avep(t^e_^c^. 

hSl^-1 


k  p 

Since,  from  Chapter  2,  "  ‘  therefore  aveptBj  •  £  A. 

yarp{Bi  •  aveptBj)  -  aveplB)  .svep(B) 

avep(3.B)  -  avep(i  f  lSi"l‘=l‘l5l‘l'=i^ 

H  k  k  It 

«  avepi  E,  E_  E.  t  c  ,t  c  .c  c  »«  o  ) 

^  pol  nti  r-1  a-1  P  p  ^  q  r  r  »  e 

^  k  k  k  k 

»  —  avepC  E,  E.  E,  t  e  t  t  o  „c  ,c  ,c  j 

‘  pel  ijel  r»i  a»6l  p  q  r  s  p  q  r  a 

I  ^  "•> 

“  c  ftC  *c  +  ^  tO  *c  »c  p  ^  iC  #c  *c  ,c 

c'p<..i  p  p  p  p  p,q»i  p  p  q  q  p,q.”i  p  q  p  q 


t  ^  I  C»c*c,  eC  ) 
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1  ^ 
avep(Bl^vep{B)  *■  c  ,e  - 

jj2  p-1  p  p' 


P-1  P'  P"‘P''1 


J 

a  Q 


•  ~?K  «C  «C  «C  +  <t.  c  *0  *0 

j^^p-1  p*  P*  P  P  p,q.».X  P  V  <3 

i  k  ^ 

^ua  varp(B)  •  i  *•“  ^  ?  i 

^  '  P  #  P  ^  P  Q  P  /» ^■*•1  p  Q  4  P  ^ 

■kep(B}  -  aVBp{(B  -  aTep{B) )»(B  -  avep[B)).(B  -  avepCfi) )) 

m  avepCBaB.B  -  BJS.avep(B)  -  B»avcp(B)  J  -  av«p(B}«B«B  +  Bw«vep(B)  »avep{B) 
+  avap{B)<,B^vep(B]  +  avep(B)  »avep(B)  ,B  -  avepfs)  .tvep{B)  .avep{B) 

«kep{B)  •>  avep{B^*B)  -  avep{B.B),arep(B)  -  avep{B.a-J-»p{B)  .B)  -  avsp{B)^vep(B.  B) 
+  2  avep(B)^vep{B)»avep(B) 

k  k  k  k  k  k 

k-" 


,  k  k  k  k  k  K 

av«p(BJB^)  •  avep(“  €.o..  S.c,  c  .) 

,3  P-i  P  P  <l«i  1  <l  r«l  r  r  s-j.  a  a  tal  t  t  u»a  u  u 


1  *‘t> 

(ir  E  c  ,c  ,c  .c  ,0  ,c  +  r"  (c  .c  ,c  ,c  ,0  ,c„  c  ,c  ,c  ,o„,e  ,c 

j^3  pml  pppppp  P».^'PP?P<14  P  P  P  <1  P  i 

PP<1PP4  P<IP?1P  SPPPPl 

^  C  •O  •C  •C  ^  C  C  nC  9C 

p  p  p  <1  «  p  p  p  a  p  q  p  p  q  p  p  q  p 

qpppqp  ppqqpp  pqpqpp 

+  c  .c  *c  .c  .c  .c  -f  c  .c  »c  «c  .c  ,c  +  c  .c  (C  .c  •<:  *c 

q  p  p  q  p  p  ?  q  q  p  p  p  q  p  q  p  p  p 

+  C  c*:  «c  .C  ,C  ,C  }  *■  £"  .(c  .C  «C  ,C  eC  *C  -f  C  .C  nC  ,C  *0  .C 

q  q  P  P  P  p'  p,qjr«l'  p  p  q  q  r*  r  p*  p“  r*  <j.*  r 

'^  c  .C  .C  ,C  aC  .C  4^  C  ,C  aC  .C  *C  .C  4  C  tC  aC  aC  aC  aC 

?  P  q  ^  r  q  p  q  p  q  r  r  p  q  q  p  r  r 

+  •C  #C  wC  -f  C  #C  *C  «€  #c  »c  +  *c  «,c  »c 

prqq  pr  prqqrp  Pq<irpr 

'T  C  •C  yC  4iiC  ^  C  •C  wC  ftC  C  #0  #0 

p  qq  r  rp  P<irqpr  P<li*<i2*  P 

^  C  «C  •€  aC  .C  *0  4'  C  aC  «C  *0  #0  *0  >f  C  .C  «C  «  ^  «C  *0 

P  P  q  q.  2"  P  r  p*  q  r*  q  p*  r*  q*  p  q*  r 


-n-tlJl  ■  ■— *^ 


&rep(S,B)*  »vep(B) 


,  k  k  k  _  k 

— (  *c  ,c  ,0  <■  ^  ,  c  ,c  .0  .c  +  I"  ,c  »c  «c  .c  +  JE"  ,c  ,c  «c  «c  VS  c 

^3  P-1  P  P  P  P  P  P  fli  q  p,<l-l  P  <l  P  q  p,4-j-  P  !1  <i  p'P-1  P  P 

1  ‘‘m 

»  S  c  *0  «C  «C  «C  aC  •»'  S  0C  aC  *0  *0  *0  4  S"  c  *0  •J  *C  *0  *0 

jj3'p-i  p  p  p  p  p  p  p.<i-i  p  p  p  p  g,  <i  p,g-i  p  p  q  q  p  p 


4-  2"  ,  C  ,C  *C  .C  aC  .C  + 


r** 

p^q-o.  P  P  4  q  q  q  P>q,r«l  p  p  q  q  r  r  p,q-a  p  q  p  q  p  p 


k  k^  k_^ 

4  S",  C  aC  aC  aC  aC  aC  4  S”  ,  C  aC  aC  aC  aC  aC  4  S"  c  aC  aC  aC  aP  aC 

p»q-i  p  q  p  4  q  q  p,q,r-i  p  q  p  q  r  r  p,q-i  p  q  q  p  p  p 

Pfq-1  p  q  q  p  q  q  p.q*™  p  q  q  p  r  r' 

Slnll&rly 

1  ^  ^1. 

avop(DaaTep(B)  .B)  •  =«  I.c  aC  aC^aC  aC  aO  4  £"  (o^aC  aC„aC„aC  aC 

v3  P-i  p  p  p  p  p  p  p<q-i  p  p  q  q  p  p 


ppppqq  ppqqqq  pqpppq 
pqqqpq  pqppqp  pqqqqp 

p,qjr«l'  pprrqa  pqrrpq  pqrrqp' 

k 


•ad 


1  * 

SV«p{B}  aaVep(BaB)  a  ^  ZLc  aC  aC  aC  aC  aC  4  T”  ( c  aC  aC  aC  aC  aC 

p  p*  p*  p*  p*  p  p,q-l'  q  q  p  p  p  ; 


ftvep(®}oftvep{B]  *avep{B)*« 


I  ^ 

■  •C  ^  ^  •»  (  ^  ^ 

^  tis^  tl  U  I)  t>  T  15  D,a»-.I'  OTiDDQQ  D9QQPP 

K’ . .  -  --  -  -  -  --  _  -- 

k._^ 

+  c^,c  ,c  ,c  ,e  ,c  )  +  5”  c  ,c  ,c  ,c  .c  ,c  1 

q  4  P  P  P  P  Pifl-jf"!  P  P  4  4  r  r 

Baao« 

1  k 

«k*p(8)  •  -?{„*s(c  .c  4iC  ,c  .c  ,e  -  3  c  ,c  ,c  .c  .c  .c  ♦  2  c  .e  .c  ,e  ,c  ,c  ) 

Ij3  P”i  PPPPPP  PPPPPP  PPPPPP 

4'  Z*  ,(c  aC  aC  „C  *0  aC  4  C  «C  aC  aC  aC  aC  4  0  aC  aC  aC  aO  aC 

Pifal'  p  p  p  p  V  4  P  P  P  4  P  4  P  P  4  P  P  4 

P>tP?P4  4PPPP4  PPP44P  PP4P4P 

P4PP4P  4PPP4P  PP44PP  P4P4PP 

4PP4PP  P44PPP  4P4PPP  44PPPP 

*  3  ®„«®„»®„«®„«®„»®„  -  3  ®„»®„«C-«c_««^»®„  *  3  ®„»®„»®_»®„»®„»®„ 

PPPP44  PP44PP  44PPPP 

*  ®«»®„»=»»®„«®^«®^  ■  ®^*®»»®„«®„«®  •®r>  *  ®„»®-»®„»®.,»®.«»®,«  ■  ®^«®^*®r,»®/.»®^»®r 

P4P4PP  4P4Pi'P  P44PPP  4PP4PP 

*  C  iiC  •C  #0  *  C  •C  *  C  **  C  yC  •&  4(C  yC 

P  V  P  P  P  ^  %  P  P  P  %  P  P  V  9  P  H  P  ^  P  P  P  P  <l 

'  C  *0  •c  •€  «C  -  C  *0  mC  *0  -  C  *0  *C  «C  mC  mC  >  C  *0  •C  «€  *0  *0 

p  p  p  q  p  4  p  p  q  p  q  p  p  p  F  q  q  p  p  p  q  p  p  q 


42c  *0  «C  4  2  C  «C  •€  «C  «C  •C  4  2  C  «C  .c  *0  •€  *0  ) 

q  q  p  p  p  p  p  p  p  p  q  q  p  p  q  q  p  p' 

^it 

p,q,r«l'  ppqqrr  PPqrqr  ppqrrq 


4  C  *0  *0  •C^«c_*c^  4  C_,C^#C^*C  #0  *C  4  C  *0  *0  *0  .C  4  C  *0  tC  .C  #C 

pqpqrr  pqqprr  prqqpr  prqqrp 

4  c  «C  «C  «C  aC  «C  4  C  mC  •C^^C  4  C  «C  *0  .C  «C  4  C  *0  aC  aC  aC  .C 

pqqrpr  pqqrrp  pqrqpr  p  q  r*q*r*p 

+  *  ^p’'^r**''p*^q*^r* "q  ^  ^'p*^r*^q*^p*^q*’'^r  '*'  ^p*^r**^q'*‘''p*^r*^q 

ppqqrr  p  q  p  q  r  r  pqqprr  P  P  q  q  r*  r 

prqqpr  P^'qqrp  PPqqrr  ppqr*q’r 


C  *C  a-'  a*'-  a'^  0^ 

?  y  (i  r  r 


(:  *e  »C  aC  sC  aC  ) 

P  y‘  q  q  r  r' 


-  ifO  - 


in  ^’C^C^C^C^C^C 

3  p,q,r-l'  p*  Q*  p*  r*  qt  r  p^  q  p  r  r*  q  p  q  q  r  p  r 


+  c  «C  «C  #c  *C  *C  -f  C  *C  «C  *0  aC  *C  +  C  *C  .C  (*C  «C  *C  **“  C_«C_*C  •€  *C  eC 

"p  pq-pq* 


n  n  n  r 


p  4  r  q  r  p' 


*  i«  ,(c  «c  ♦  c  .c  )«(c  ,c  4  c  ,c  )«(c  .c  4  c  «c  ) 

3p,q,r-l'  p*  q  <l  p'  '  q  r  r*  q'  ^  r  p  p  r' 


elO|5{B3  «  av«p{(B  -  av«p(B)  )«(3  -  Avept&l  )*(B  -  avep{Bl).(B  -  avepfBl  )' 

-  3  avep((B  -  avap(B) )a(9  -  avep(B) )l«avep{(B  -  avep{B) ).(B  ^  avepCfi) )) 


ave{(B  -  av«p(B)  )«(B  -  avap(B)  )«(S  -  avep{B)  ).(B  -  avrsTpCB)  )) 

■  avep(B«E.B«B]  •  avep{B.B«B«avep(B] )  •  avep{B.B*avep(B) •&}  >  avsp(B.avapCB} «B<B) 

-  avep{avep(B] •B.B«B)  4  avep(B.B.avep(B) >avep(B) )  4  avep(B.avepiB) .BaavepCB) ) 

4  av8p{B«avep(B)»avep(B)4.S}  4  aveptavep(B)  .B^.aveptB) ) 

4  avep(avep{B}  •B.avaptB}  .B]  >  avep(avep(B)*avep(B]*B,>B] 

-  3  avap(B)»avep(B).av«p(B).avep(B)) . 


After  conaldarable  algebraic  aaclpulatioo  the  fourth  cunulant  of 
B  averaged  over  randoaizatlODe  within  blocks  can  he  written  as  follows: 

•lop(B)  -  ^p,i"r,,^KCp.C^4C^.Cp)  .( VC^4C^.C^).(c^.C^4VC^).{o,.Cp4Cp.cJ 

♦  ('p-^q  V'q^'^V's  ^  's’^r^ 

4  (c  ,C  4  C  ,C  )«(c  ,C  4  C  ,C  )»(c  ,C  4  C  .0  ),(c  ,C_  4  C  *0  ) 1 

P  4  4  P  r  B  8  r'  '  s  p  p  s'  '  q  r  r  q''* 

•  ^  +  c^,c„).(c„.c^  4  c„.c  )«(c^.c„  4  c„.e^))  4,1 

p^q«i'  rq  qp  pq  qp  pq  qp  pq  q  p 


for  reference  we  repeat  here  the  results  for  the  lower  cuisuitots,  cojaely 


3kep{B) 

1 

e  ^  2  _(c  *0  +  C  *0  )*(C  mC  +  C  #0  )#(c  +  C  *0  ) 

3  p,q,rel'  p  q  q  '  q  r  r  q^  '  r  p  p  r^' 

4. a 

varp{E3 

X  j'l*  f  \ 

“  ~2  p7q-l'''p*"q*"p*"q  ''p*''q*''q*‘'p' 

4.3 

avep  (b  1 

1  j;  1  . 

»-r  ^-C*C  »  -rA 

k  p»l  p  p  k 

4.4 

-  - 


Chapter  5 

DfRIVATIOS  or  OTS  CU>«nAaT3  or  »  AVXRAOID  OVBH  B0TA?IC«3 


1 ,  Soaie  Introductory  jtoagrjt* 

Wo  ere  Interested  In  aesre^ng  the  randtxlBetlori  (perautetioii) 
cuaulusts  of  B  over  rotatl<n>s  of  the  coordinate  systen  In  vhlcb  the 
vectors  defining  B  ere  expressed.  The  derivetlac  given  below  is  for 


the 

2>dlaenalonal  case. 

Let  the 

vector  Cj  be  denoted  by 

{•i.  »>^) 

for 

1  <  1  <  X 

vbare  the 

and 

bj^  are  scalars  for  any  1 

.  BscaU 

now 

that  If  B 

la  looked  at  In  a  coordinate  systea  such  that 

A  Is  of 

tbs  fora 


1  0 
C  J 


then 


^  2 
Z  af 
1-1  1 


Jc  2 
Z  K 
1-1  1 


and 


X 

Z  a  b 
1=1  *i”l 


Consider  nov  an  arbitrary  rotation  of  axes  In  the  (a,  b)  plane 
such  that  the  nev  coordinates  are  denoted  by  (a',b'}.  Then  a  point 
(aj:,  b^)  In  the  old  coordinate  aystea  will  now  be  given  by  (a^,  b^) 
In  the  nev  coordinate  aystea  vbero 

a^  =  SOS  S  =  »ln  S 

b^  »  a^  sin  e  ♦  b^  cos  9 

and  S  Is  ths  angle  betveen  the  old  and  the  nev  axes. 


2.  Derivation  of  the  Cwaulcnts  of  B  Averaged  over  Rotations 


Frco!  Chapter  4 


I  k 


aveyCE} 

*ve  f  r  2, 

( rotat-iona )  ^  1  ' 

-  1;^  r.  ;  (a!®  =! 

^2* 


iZ 


,  2*  k 

^  i 

0  l*u. 


^  r  ‘- 

^  ‘ 
&>  ^ . 

A 

5»l  i  1 

k 

k  2 

i5.i< 

Si®!®! 

-  b^  Bln  S)^de 

o 

1  +  b- 

eln^e)de 

.  S  ^  k  2*  k  p2«  2 

■  3-{  E  /  co»‘'8d9  -  2  E  a'b'  /  cos  9  sin  e  dS  +  t  it! '  J  aln  9  de] 

1«1  ^  0  1-1  ‘  ^  0  1-1  0 

ki  X  li  ff  ^ 

•  Ij'  £  aJ^  /(i  *■  icoa  2d)iie  •*■  |^  C  /cos  ed(coB  S)  +  —  £  b^'^/  (•J  -  ^cos  2e)dfl 

i-l  0  "i-1  0  1*1  0 


k  „2* 

£  a 

’'i-1  *  0 

k  ,  k 


-  i  X  a!‘=  f  i  £  b  • 
1-1  ^  1-1  ^ 

k  Q  O 

.  A  E 

’  iJi'*!  l' 

illao 


,2 


r  1.  n 

(rot.)  i-l*!**!  *  ° 


and 


*v«.  T  2  ,  f  ,  2  P. 

(rot.)  l^**!  *  4  l5l(»l  + 


ave  a\ep(B) 
(rot) 


\ 


.  »  O  -J  / 

*  t  i  3“  +  b" )  / 
•i-il  i  1  7 


ave  varp(B)  avc  1  ,  j—  ,  ,  „  ,11 

(rot)  “  (rot.)  ^  P(?*l  P*”<1  P 


C  «C  aO  «c 

p  q  p  q 


(a  ,  b  ),(a  ,  b  ),{a  ,  b  ).(s  ,  b  ) 

p  p  q  q  p  p  q  q 


aa  ab  aa  ab 

,  p  q  p  q  >  ,  q  p  q  . 

“  (  -I 

^  b  a  b  b  'b  a  b  b  ^ 

pq  pq  pq  pq 

Thl*  product  cao  be  written  as  a  <i  x  >.  b  x  P  ari-ay 


aaaa  aaab  abaa  a  lab 

pqpq  pqpq  I'qpq  pqpq 

1 

,  I 

aaba  aabb  ■'  ^abba  abbb  '' 

pqpq  pqpq  pqpq  pqpq 


baaa  baab 

pqpq  I'  q  p  q 


b  b  u  a  b  b  a  b 

p  q  i>  q  pqpq 


b  a  b  a  b  a  b  b  ^  'tuba  u  b  b  b 

pqpq  p  q  j'  q.  ■■  q  p  q  pqpq 

Tblfi  product  of  two  dyavia  contains  eitiuinta  of  six  CLtstlnct  tyiWE  ; 


(1)  a^'a*"  ;  (P)  a"a  b  ;  (3)  a  b  a  b  j  (U)  aTj'  j 

pq'  pqq'  pqpq'  pq' 


(^)  a  b  b~  i  (6)  b^b^ 


j,  the  product  ol'  d^ads  in  the  variance  of  B,  can  aino  be 
writtea  as  a  exnxext:  army  anu  the  distinct  eiementa  appearing 
in  Buch  an  array  are  Ideoti-.-al  with  those  above.  We  /ill  now  coneldcr 


the  e 'feet  on  ti*  elements  (1)  to  (6)  when  the  coordinate  syatem 
(a,  b)  ia  rotated  thix'uah  an  arbitrary  angle  e  and  then  the  average 
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in  over  O  <  £•  <  2, 


k  2* 


P  «  2s 


Since 


p,q"i  0 


k  2« 


i  y  •'  k^sli".  e)‘'(a^C08  0  -  'B’sln  £>)''dS 


2^  E-  /  (^;-co.-e  -  2. 

P,q-I  G  P 


^pb^coe  e  Bln  a  r  V'  ljin^’a) 
9  Bin  9  + 

D.C^I  o  P  1  P  Q  a  ' 


■"e  *in  e 


*  '**p*’p*»a^PP*  S  •in^S  +  b'S'^Bln**! 


p  a 


s)de 


1  f*  If  I  2i!  , 

2,  /  C38  9  da  .  „  /  d0  - 


Si 


^  I  cos ^9  alu%  AS  m  i  4-d 

0  O 

1  ^  gjj 

/  CCS-  9  Bin  e  as  .  i-  /  COB  0  Bin3. 


B  AS  m  0^ 


one  obtains 


•ve  2  2 


{rot)p,|:;i'‘p»P  -  ^  ^  ^ 

■  -  n  lb'2  .  i  ^  /..2  .  ..2.2 

-  F  P-i  p  o  p»l'-p  '  "p  / 


(a  +  b  3*^ 

p»i'  p  p' 


b«c*.iiae  a  '  +  "b  ’  ' 

P  ? 


2 

a  t  b 


and  Sa. 


P 


Eb 


T>' 

P 


(2) 


RVO  j-..  2  , 

(rot)  p/-i-0.*p‘q  q 


(3)  r  .  b  a  b 

(rot)  p  p  q,  q 


^  (rot)  p,q«l  p  St  “  "5  p«l^  p  ^p^ 


-  r  /  2  v2v 

5  p»l'  p  p' 
2,2 


(5)  /“Ifs  ,a  b 

(rort)  p,q«a  p  p  q 


C(»aldsr  ncm  the  third  cuouleiit  m  given  by  Tonoul*  ^  .2 

8 

Bkap(B)  •  £"  ,(c  +  C  «C  )(c  sC  4-  C  «C  )(c  aC  4  c  •«  ) 

v3  P  q  q  V  q  r  r  q  r  P  p*  r' 


*  ^P.qpr-l(V%-S-"r*%*<=r  "  V"q*=p“=r* V't  ^  V\'V 
♦  “p-'q-^'r'^p-'r'^q  ^  ‘^p*=q-‘=q*°r«‘=p**^r  ^ 


pqrqpr  pqrqrp' 


By  a  slallar  procedure  me  vae  used  for  the  averana  over  rotatloas  of 
the  variance  of  B  It  can  be  ehown  that  all  the  above  teme  is  the 
■kevneas  of  B  cootaln  only  eleaents  of  the  following  typeas 


-  46  - 


(i) 

r»  2  2  2 

p,q,r-?L 

{0 

r  .aW 

per 

{2) 

V"  2  2. 

JS  Ig,  •  * 

n) 

b 

.a  b  a  b  1 
P/4;r«l.  p  P  4  q. 

(3) 

E" 

P,<i,rmlp<ir 

(8) 

k  ^ 

^  ,a  b  a  b  t 
p,q,r«l  p  p  4  % 

w 

P/I.r-dVl'*  <l^r\ 

(9) 

r"  -w  v2v2 

(5) 

P»4jr^*p“4\^r 

(10) 

r  ,bW 

p,(l,r«a  par 

Zcatead  of  doing  the  rotatloo  og  thase  eleaenta  directly,  m  mb  done 
for  the  first  tvo  cuBulaot*,  an  altarMtlve  procedure  vlXl  be  used 
vblcb  ellxlaates  aoaie  of  tbe  drudgery  of  the  algebra.  Tbia  proceoure 
proree  to  be  very  aetlafactory  for  the  av'ereglng  of  the  fourth  ciBsolaDt. 

'Bm  •Itemative  procedure  aakea  um  of  the  feet  that  after  averag- 
log  over  rotations  the  eleaenta  are  given  by  lioear  coiAlnatlcna  of 
the  Invartants.  By  finding  all  the  Invarlaota  In  advance,  or.j  can 
•aelly  find  the  coefficlenta  vhleh  expraae  each  of  the  elcaents  as  such 
linear  cosblnations. 

91  tics 

a  •  a'  cos  9  •  b'  sin  d  and 

h  -a'  Bin  e  ♦  b'  cca  9 

then 

a=  a  a  cos  a  V  0  sin  9  and 

b'  «  -a  sin  9  ♦  b  cos  9, 

so  am  InfinitesisiaX  chaunge  dd  in  9  vlU  change  a’  by  an  aiaount 
b'l^.  Banc«!  an  InfinitesiaHl  ebauge  A9  la  9  vill  change 


E  a  by  an  amouiit  fx'v?  £'  ,a  b  n  o'  .  A  two-way  table  can 

P,q,!--A  p  q  r  •'  P,q,n.i.  p  p  q  r 

b«  com-; liuetbd  to  shov  the  effect  of  an  inf InltcaliripJ-  change  In  6  on 
S“ch  :if  ths  ten  sl-siasrits  t-o  the  skevness. 


Tsole  ^.1 

CoefficlentB  of  the  Linear  Combinations  of  the  Kleroenta  of 
ekep  B  aft.er  an  Infinitesimal  Cliange  Id  $ 


7‘h«  mnnher?  in  TS-hls  5.1  refer  to  tho  coeffloi^rtito  in  the  iir,»aj- 
cojiblnatl  one  obtained  by  itaXlng  an  infinlteslmel  change  in  f<,  Ftr 


K,,  HO 

example,  a  change  fiO  In  h  changes  S  ,a  b  b  “o”  by 

p,q,r«i  p  p  q  r 


k  o  o  '-i 

( -  r”  -  a^b“b' 

■  P>a,r«=i  p  q  r 


k 

t  r"  «  b  ! 
P,q,r»a.  p  p 


V.  .2 

q  q  r 


k 

r" 


p,q,r»l"p“q“r' 


To  find 


i*6  - 


«Jj  the  invariaots  under  rotetlon  ane  or4ly  nstjds  to  find  thoee  ctwo- 
bluBtions  of  (  i  )  throogt'.  (10)  vhich  are  unaltered  by  an  Infinitaoiml 
chants  In  -i*  *  itiupe  a/  e  i 


(a)  -  E"  +  2  a^a  B  a  b  +  ■!-2ababb^  + 

p,(l,r«J2  p  q,  r  pqqrr  pqr  pqr  ppqqr  pqr' 

m  E"  ,  [a^(,.'’6‘‘  +  2  a  b  a  h  +  4-  b‘^(a‘"a‘^  +  2  a  b  B  b  +  bSs^) 

P,il,r-^  p'qr  q  ^  ^  qr  P  qr  (iqrr  qr' 


E"  (a^  +  b^)(a  a  +  b  b 
;,r-l‘  p  p  q  r  q  r' 


P,q.-r' 


«  Z  (a^  +  b^)  E"  (a  a  +  b  b  )^  -  £X",(a®  +  b^)(a  a  +  b  b 

p»i'  y  p'Piq-1'  p  q  p  q'  p»q"0-  ?  p^'  p  q  p  q' 

But 

r  ,(a  a  +  b  b  )^  -  [(  E  a*)^.  E  a^)  +  2{(  E  a  b  f  -  isV)  +  [(Eb^)-  -  Eb^] 

p,qmV  p  q  p  q'  "p*l  p'  p«l  ‘'■p-1  p  p'  p  p'  p" 

.  2  -  I  (."  . 

P-1'  P  P' 


axiil 


E"  -{a^  ♦  b^)(a  a  *  b  b  )‘ 
P,qal'  P  p"'  p  q  P  q' 


^  L 

Eit  f  H  ,v. 

B  (a  a  > 
p  ^ 


2b, \  a  b  +  +  2a  b^a  b  +  b**a^) 

ppqq  p.?q  ppq  ppqq  ?q 


r,(a^  -  a^  4  2a\^'  4  a\"  -  a^**  4  aV  -  aV  -  2aS'-  4  -  bf ) 

p-1'  p  p  PP  PP  PP  PP  PP  PP  P  P 


=  .  2  .2.2  2  Jii3 

E  (ft  +  b  )  -  E(ft  ■♦•  b  I’i 

pal'  p  p-"  '  p  p" 


llauB 


(b) 


.  4,9  - 


(») 


SlBlltrXy 

k 


1*  -  2  r.  +  b®)‘’  -  2  E  (b^  +  -b^)^  +  2  E  (»^  4  V 

p-i-  p  p'  p«a.'  p  p'  p«i'  p  1 


4  .kE(.^.b2)2.2E(a'^.b2)3. 

r«i  p  p  p-i  p  p 


2^3 


.  8  -  6  E  (»2  +  b^)S  ^  2  E  (a®  ♦  b^) 
P«1  P  P  P^  P  P 

and 


(c)  -  E"  (a^a  b  a  b  -  a^aS>^  +  a^a  b  b^  -  a^S®) 

'  ‘  P/Ij^'*!'  pqifr  par  P44r  pqr' 

-  ^  ,(abab  -  a\®)(B^  +  b®)  -  -2  +  E  (a^  4  b®)^. 

P><l4r«a'  P  P  4  a  P  r'  r  r'  p-l'  p  ?■' 

Tbaae  tbrea  iDvarlaots  (a),  (b)  and  (c)  are  such  that  00  other 
linear  coablnstlAD  of  (1)  through  (10)  can  be  found  Independent  of 
thea.  Purthencore,  these  three  Invariants  are  Independent  of  each 
other  and  are  «lLL  expressible  In  the  form 


a  4  p  E  (8*=  +  b^)^  7  4  b^>- 

P*»l  P  P  'p-l'  p  p' 

vbere  Cl,  p  and  y  are  constants.  Consequently 

1,  (a^  4  b^)®  and  E  (a^  4  b®)^ 

'  P'1'  P  P'  P*l'  P  P 


k2v3 


waa;  a.  V  -  wai  V  C  li  A  « 1  i 


iVaxl^ivo  auiu  exeSKni'^if  Or  Xji€  XOl^Xl  CUMU* 

lanta  of  B  ( vhen  averaged  over  all  rotations  of  the  coordlne.te  system) 


2  2  2 

must  be  expressible  as  a  llnetir  coohinatlon  of  1,  £,(a  4  b^)  and 

E^(a^  4  b^)3. 

P-l  P  P 
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By  choosing  suae  einple  sets  of  vectors  one  can  easiJy  compute 
the  values  of  (1)  through  (lO).  This  laethod  vili  nov  be  oiitll.-jel. 
let  US  choose  the  fnliowlng  three  sets 


vW-«  ■  -k/A  t>  • 


(a) 


O) 

ji/Js  0  I 
i/y?  0 
0 


1/y^ 


0  i 

o/ 


for  the  ast  a 
k 

P-1  P 

for  the  Mt  P 


'  Th  «  ' 


2, 


Ma"  t  b2)3 

P-1'  P  p' 


r  /  2  .2»2 

£,  (s  eh) 
P-1  P  P 


1, 


P-1  p  p 


for  the  set  7 


r  I  2  .2.2 

E,  (a  +  b  ) 
P-1  P  P 


2 

3' 


k  p 
(a 

P-1  P 


Any  third  cunulant  eleoent  averaged  over  rotaiiona  la  of  the  fom 
-  t.yt(a2rbj)^t  C£(a2tb2)3 
idiere  L  la  an  arbitrary  expreaslon  of  sixth  degree  lo  a  and  b. 


“o»  Po' 


7^  be  the  values  of  an  arbitrary  quantity  L  averaged 


over  rotatlKis  for  the  sets  a,  ^  and  7  respectively. 


a  •  X  r  2y  +  £z 
o 


X  *  y  +  s* 


-  *  +  3  y  +9« 


Then 


and,  solving  these  equations 


aveC,J 
( rot . ) 


-  *4®  +?v)  +  {-j^  tflp  -  t  t! 

't  o  o  i:  c'  'ho  o  h  'oTwl'  p 


2.2 

1? 


+  {|-a  -6p  +  2.j,)E  (a,^  + 

'  2  o  o  2  'o  p-1  p  p 


A  t*bl*  of  a^,  and  7^  lor  (1)  through  (lO)  Is  •bown  below 

AS  Table  5.2.  The  values  of  the  e<5efficlents  x,  y  sand  t  ixi  the 

expression  for  sretl}  osy  then  be  wltten  down,  using  this  table, 
(rot.) 

nMse  values  are  given  in  Tftibla  5-3. 


Table  ?.2 

Velues  of  the  eleinonte  in  skep  B  averaged  csver  9  for  the  rectors 

a,  fi  end  7 


'table  ^.3 

CoefflclentB  in  the  linear  coBitinatlonn  expressing  the  eleaents  of 
ave  step  E  in  terac  of  three  convenient  Invaxiaats, 

(rot)  _ _ _ _ 


- - - A  -  X'  " 

■  >>  .  .p' 

- - jrTV 

r(«"  -I 

L  UP  :p' . 

(i) 

1 

-9/8 

5/8 

( ‘  - 

0 

0 

0 

(3) 

1 

-5/8 

1/3 

(h) 

0 

-1/0 

1/8 

(5) 

0 

0 

0 

(6) 

1 

-5/8 

)/8 

(7) 

0 

b  ' 

b 

0 

-i/0 

1/8 

(9) 

0 

0 

0 

Llol__ 

i 

-v/6 

_ 2/8 _ 

The  fourth  oui/iulant  of  B  contains  product  of  '^vmda  of  the  1‘orir 

c  .c  ,0  ,c  ,c  .c  ,c  ,c  and  of  the  form  c  ,c  ,c  ,c  ,c  »c  ,c  ,c  , 

P  Ql  P  r  q  u  r  B  P  11  P  q  p  a  P 

To  find  the  average  over  rotation*  of  elop  B  we  first  need  to  find 

a  moxlmfU.  set  of  Independent  Invariants.  It  will  then  he  possible  to 

find  that  linear  combination  of  these  invariants  which  results,  after 

rotation,  from  each  of  the  distinct  eisaentu  in  c  ,c  .  . . .  «c  o.id 

p  q  8 

c  .c  e  ...,o  .  The  distinct  eleinenta  arc 
P  <1  q 


1) 

i  p  ?  2  2 

p,q,r,3v.i  p  q  r  8 

16) 

k._ 

p,q«i 

it  1 
a  a 

p  q 

2) 

d  ?  2  , 

a  a  <31  a  V 
p  q  r  s  » 

17) 

it  3. 
a  8*1) 

p  q  q 

3) 

p  q  r  8 

Ifi) 

V  2,  £ 
a  a  b 

p  q  q 

2  2  »  V 

a  ft  e  b  a  b 
p  q  r  r  0  a 

19) 

**  1 3 

a  a  b 

p  q  q 

5) 

' '  p  2 

a  a  'a  b  b 
p  q  r  r  0 

20) 

‘‘v'i 
a  b 

p  q 

6) 

p  q  r  s 

31) 

a\ 

p  p  q  q 

7) 

a‘‘e  b  a  h  a  k 
p  q  q  r  r  3  fa 

.:>) 

aV 

p  p  ■!  q 

3) 

a“8  b  a  b 
p  q  q  r  r  a 

•5  "5 

r.-^b  a 

p  p  q  q 

9) 

o  J  j 

»“&  b  n'T)'" 

;•  q  q  r  B 

2J.) 

a^b  b^ 

p  1  q 

10) 

p  q  r  s 

25) 

2.  2  2.  ? 

a.  b  a  0 

p  p  ^ 

li) 

a  n  a  t>  a  b  a  b 

p  P  q  q  r  b  3  s 

^6) 

a'^bi^a  b'^ 

p  p  q  q 

■! 

*  ; 

,  ,  ,  .  <1 
d.  u  a  L?  a  p  B 

P  C'  q  q  r  r  5 

27) 

o  o  li 

a~b  0  ■ 

p  p  q 

lU) 

a  b  a  b  b'^b 
p  1  q  q  " 

a  t'  b‘'b‘  o' 
f>  p  q  r  s 

rf) 

a  fc  a  b' 

P  p  q  q 

a 

p  P  q 

15) 

b  ' 

p  q  i'  i 

iCj  ) 

4  4 

b  b 

P  a 
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TrOi\  .art  ji'  Table  I  nv-... .  5/ 1  ng  (1)  tlircugn  {-I*"')  the 

ro.j  owing  i-et  of  Independeet  Invariants  can  be  deduced; 

(a)  «  (1)  +  ^(4)  ■*  i^(6)  t-  U(il)  +  4('.3)  + 

(B)  ..  (1)  +  )*(3)  ^  £(6)  ,  +  (15) 

(C)  -  (6)  -  ;?(0)  +  (11) 

By  rearrmnaenent  and  soma  algebra  (K),  (B)  end  (C)  can  be  vrlttea  aa 


8  I  (a^  b^)->  8  i  e  -  1.  e  h^)^* 

v'l  ?  p'  p-1'  p  p'  p-1'  p  p' 


B  -  16  -  6j^(a2  e  b^)**  .  l6j,(aj  e  bj)3  -  2bJ^(,2  , 


C 


r  /  2  vl\2,2 

P-1^  P  ' 


is.  ia 

1  -  e  *  I  r  ,(a  b  -  b  a 

p«^'  n  n*  O  Tx  #a«i  '  T\  A  T\ 


p»i'  p  p’  2  p,q«l  p  q  p  q' 

Tram  the  table  Involving  thrcsigh  (30)  the  following  set  o' 

independent  Invvtrlants  can  be  deduced 

(A')  »  W")  +  M-il)  +  6(25)  +  1(28)  +  30 

(B')  -  !,i6)  +  i(l6)  *  2(20)  e  e  *^(27)  +  (30) 

(CO  -  (20)  -  ‘*(23)  *  3(25) 

Theee  can  also  Oe  written  ae 


(a‘  )  -  E"  ,  (a  e  +  b  b  ) 

p,q«i'  p  q  p  q' 

(B')  -  [  V  u'")"'!'-  -  i;(a"  b“)^ 

'  “p*l’  p  p-  '■  p  p- 

1  j. 

L  1  ^  .J  ^  a  Q  —  y  a  j 

p,q*.I'  p  I  r  Q 


froi7.  a  careful  study  of  these  invariantB  it  soon  bccoines  apparent 
frir-t  a  ccr.venlerit  siaxlnyal  eet  of  Indepf.ndent  invarleuite,  as  linear 
coniblnatlona  of  which  we  may  express  the  average  over  rotations  of 
the  fourth  cowulant  of  B,  are  1, 


p  p'  ’  pH'  p  p 


r  I  2  .2.4 

p«l'  p  p  ' 


t  ^  .(fc  +  b  b  )* 

P»1  P  P  P><1“1  P  9  P  9 

Q%e  tivemge  veuues  of  the  terras  (16)  throutCh  (30)  for  rail  rotratloos 

of  the  coordlprate  syitea  only  iDYolve  the  llnerar  con^lnratlona  of  thrM 

IswrlantB  (A'),  (S'),  (C).  A  trable  of  thb  coefficient*  of  {A’)^ 

(B'),  (C)  for  eecb  of  the  fifteen  le  now  given  rae  Teblo  5.5, 


r  {,  b  -bra  )“ 
P>9"41  ?  9  P  9 


Table 

CoefflcientB  of  the  llnerar  cosfctnratioos  of  (A'),  (B'),  (C)  for  the 
eleaests  (I)  through  (l^)  averaged  over  rotation*.  (Bach  entry 
naut  be  divided  by  12&) 


(A') 

(B-) 

(C) 

(A')" 

(B-) 

(O’) 

(16) 

20 

.45 

-12 

(23) 

0 

3 

-8 

(i7) 

0 

0 

(24) 

0 

0 

0 

U6) 

-4 

9 

-4 

(25) 

4 

-1 

4 

(19) 

0 

0 

0 

(26) 

0 

0 

0 

(£0) 

-12 

15 

20 

(27) 

-4 

9 

-4 

(21) 

8 

-3 

0 

(28) 

8 

-3 

0 

° 

0 

0 

(29) 

0 

0 

0 

■ 

_ 

(30) 

20 

15 

~  56  - 


The  R/erage  value*  of  the  terms  (i)  thrcijgh  (15)  for  esj,l 
rc  tatlons  of  the  ooortilnate  syateii  involve  edl  of  the  Bevea  Invari- 
ant*  aad  are  given  In  Table  5.6.  Both  Tablee  5,?  siid  5*6  were 
calculated  by  direct  rotation  method*  as  given  earlier  in  this 
section  for  the  necotid  cunulant  and  alno  by  the  neth^d  outlined 
for  the  third  cumilant. 


Table  ?.6 

Coefficient*  of  the  linear  coohlnatiocis  expressing  the  element*  (1) 
through  (1?)  in  terms  of  a  set  of  invturiaats.  (lach  entry  must 
be  divided  by  126.) 


D 

ISIS 

lUjUl 

i:(a^+h^)^ 
^  P  P' 

Ak}>  >b  a  f 

1  P  V 

(1) 

1^8 

-288 

320 

-150 

1*5 

60 

-36 

0 

0 

0 

0 

0 

0 

0 

128 

.192 

~k2 

27 

-12 

-12 

(V 

0 

-16 

-18 

3 

12 

-4 

(5; 

0 

0 

0 

0 

0 

0 

(6 

128 

-160 

a 

-22 

13 

-4 

26 

0 

0 

0 

0 

0 

0 

C 

Wm 

0 

-16 

32 

-lit 

5 

4 

-12 

9) 

0 

0 

0 

0 

0 

0 

0 

(10 

126 

-192 

128 

-42 

27 

-12 

-12 

(U. 

0 

0 

0 

-6 

-3 

12 

12 

(12 

0 

0 

0 

0 

0 

0 

0 

13^ 

0 

-16 

32 

-18 

3 

12 

-4 

14 

0 

0 

0 

0 

6 

0 

0 

15 

128 

-288 

320 

-150 

**5 

60 

-36 

By  codblnlng  the  elements  (l)  throu^  {1^}  with  elements  (I6) 
through  (30)  according  to  formula  we  arrive  at  Table  5.7  which 

has  ^1  entries.  Some  one  of  these  21  terms  la  eijui valent  to  each 

Q 

of  the  256  terms  which  arise  in  the  2°  array  defining  the  fottrth 
cuaiulant  of  B  a/emged  over  rotations  of  the  coordinate  system. 


,4  a 
^  t 


•w 

P  4  ^ 

As," 


a\b3 

•M' 

P  P  4  4 

A  A' 

P  P  4  4 

•^pVJ 

P  P  4 
P  P  4  4 

,A"A" 

P  P  4  4 

•S,"«  b^ 

P  P  4  4 

P  P  4 

P  P  4  4 
2.2  .1 
ft  O  »  B' 

P  P  %  q 

aS)S)^ 

P  P  4 

a  b^a  b^ 
P  P  4  4 

a  bS,** 

P  P  « 

Vq 


W  t"  ,aSW 

p,*,r,i««l  p  4  r  e 
2  ?  2 
*P*4*r* «'*« 

aV.V 

P  4  r  e 
2  2 

a"."*  b  b" 
p  q  r  r  » 

a"aS.^" 
p  q  r  s 

2  2 

p  4  ^  r  8  • 

2 

4  4  b  4  b  e  b 

p  q  4  p  r  « 

•"»  b  bS," 

p  4  q  r  » 

aSi"»  b  a  b 
p  q  r  r  (  > 

."a"a"b" 
p  q  r  B 

aS)"a  b  b" 
p  q  r  r  8 

a"bS.S." 

p  q  r  s 

a  b  a  b  a  b  a 
p  p  q  q  r  r  * 

9 

a  D  a  D  a  3  b 
p  p  q  q  r  r  * 

a  b  a  b  bS," 
p  p  q  q  r  9 

2  2 

ft  b  b  a  b  b^ 
p  p  q  r  r  9 

a  b  b^S:" 
p  p  q  r  8 

p  q  r  8 


Chftpter  6 

JtORMAL  TILSOHf  FOR  TWO  T5(iCA'»SKH’ni 


u . X  iat-rciCLuctlon ;  t*  In  Oae  Plaenslon 

A  BtaadArd  laEthod  for  testing  whether  the  loeen  of  a  najuple  has 
a  certain  specified  value  when  the  population  variance  is  unhnovn  Is 
to  Bake  use  of  the  t  distribution.  Suppose  that  x^,  x,^, 
la  a  sample  of  k  Independent  observations  from  a  nopoal  distribution 
with  unknown  aaan  p  and  itandard  deviation  a .  To  teat  the  hypo¬ 
thesis  that  tbs  aean  of  the  population  has  some  speclf.,.ed  vsil’ie ,  say 
u^,  one  forma  the  ratio 

Ll-ffi.  6.1 

m/Jk 


k  /  '  k  ^ 

where  IT  »  and  a  ■  x)® 

This  ratio  has  Student's  distribution  with  (k-l)  degrees  of  freedom. 

In  the  asjiK*  way,  a  test  for  the  difreronce  between  two  population 
means  makes  use  of  the  assumption  that  the  populations  are  independently 
nonrcilly  distributed  with  unknown  means  p,  and  u-.  Tf  the  variances 
are  also  unknown  but  can  be  aBEumed  equal  then  a  statistic  for  testing 
vR;»ther  the  ac?eri3  are  equal  uakeu  use  of  a  pooled  estimate  of  the 


pjopuJ  at  I  on  vorl.aiice  dr^flned  as 


2 

s 


-  1  1  --  .“^rir 

'ij'-a 
X  o. 


1  ^ 
-  / 


where 


H 

(i)  o;-  ^  (x.  -  3c  i  ^  i,  ^ 

i  .V  j-i  J»J  ij  1 


(S6«f^e4f  4Wiss«4si'«*-  ,  s.-i-'r-;;-.. 


“li'  "i:-' 

SfUJ::.-'  . 


fii'H  OilO 


(^-1.0 


ui  thtt  poptiiationb 


i  Kj  j-i  i; 

The  statistic  used  i» 


6.2 


The  ratio  Is  distributed  as  t  vitb  (k^  *  ^2  ~  degrees  of  freedon 
and  the  test  la  made  In  the  same  wajr  as  the  single  ssaiple  test  given 
above 


If  one  Is  not  able  to  make  the  assunqptlon  of  equal  populatloR 

variances  then  the  problem  Is  more  ccui;>lex.  The  separate  estimates 
2  2 

Sj^  W)d  Sg  for  the  population  variances  are  used  In  the  ratio 

6.3 

This  ratio  Is  not  distributed  as  t  with  (k^  *  kg  -  2)  degrees  of 

freedom  as  vas  the  case  vlth  equal  variances  because  in  general  the 

denominator  la  not  prop<)r,;lonea  to  «  X  with  (k,  +  k,-,  -  2)  degrees 

1 

of  freedom.  Mo«ev^  ,  this  distribution  has  been  studied  by  Welch 

(16)  and  algnificancB  points  have  been  tabled.  KUrtheraore,  Welch 

[17]  showed  tfiat  the  tjereentege  points  can  be  approximated  by  the 
percentage  points  of  the  t-  ■!,  strlbutloa  with  f  degrees  of  freedom 


where 


k  ,  ■  h  '■ 

1)  -f  ■  j)  --J^l 


i  *;  X  i  .  iT. 

TTiene  niake  U-se  of  the  ab6un;ptioi)  that  the  uiiJerlying  popu¬ 

lations  are  tioraioliy  distributed  but  as  mentioned  in  Chapter  1,  they 
are  remarkably  insenal'kiva  to  deviations  from  nornia-llty . 

2 

6.2  Introduction:  Hotelllng'B  T 

Proceeding  now  to  the  multivariate  problem,  we  assume 


lU 


U'  *21' 


Ifc' 


are  pairs  of  observations  drawn  Independently  from  a  bivariate  normal 

2  2 

dletrlbutloo  with  paraaneters  p, ,  c,‘,  Og  and  p.  The  sas^rle 

means  (x,  ,  x,  )  are  distributed  accord,  jg  to  a  bivariate  normal  dls- 

Jl  d 

trlbutlon  with  parameters  u^,  c^/k  and  o.  These 

2  2 

sample  means  are  Independent  of  s^,  and  r,  the  estloates 


of 


Og,  and  0  recpeottvely  which  are  defined  as 


,  k 

k-J  ’«!  ij  I' 


)(xo.  -  k  ) 

■  “  ■  --ji _ _ _ 


i,  2. 


6.3 


.  k 


The  quantity 


!J:,  )' 


(*i  -  M,  )  (x,^  -  p  ,)  (x,j  -  p,^)  q 


i.-P 


a,  /K 


d^Jk  ] 


It;  dlatributvd  at  X  with 
(.•“l  and  p  by  the  estiiaates 


degrees  of  freedom.  Replacing  a^, 


6,4 


6.4  i replaced  by 


and  r  respectively,  expression 


6.5 


i _ ^  ^ 


Wheij  *p<)j  ^  i  .1  *'?  iridependent  aamplcR  from  ft  fciv!!..rlat,e 

*2 

corml,  then  T  is  distributed  as 

"  fi-i#  b,K-S 

To  test  the  h^TpotheBlft  *t  the  0$  level  that  aunti  have 

speell’led  raluea  and  say,  against  the  alternative  that 

^  and  Jig  /{  one  cenputea  T  vlth  and  replaced 

■fay  and  and  ccaaparcie  the  result  with  the  value  of  6.6  using 

the  table  of  f  with  2  and  (5c'2)  degrees  of  freedom  at  the  Qji 
point , 

The  two  sairplft  prob'lew  In  .sore  than  one  dlrenuloB  -follows  a 
eladlar  procedure.  This  will  be  iUuatrated  by  tenting  that  the  dlf^ 
ference  between  the  means  of  two  blvaxlats  normal  dlstrlhutloaa  ml^t 
reasonably  bo  tero.  het  and  (y^^,  y^^),  1  <  1  <  1^, 


ujsd  1  <  J  <  «.£,»  he  independent  aaaples  from  blva- late  r.onnsd  popula- 

2  2 

tlons  with  pareunetera  p  a  ,  a  ,  p  and  p  ,  a  , 

*1  *2  *1  *2  *1*2  yi  yi 

2  2 

Oy^,  Oy^,  *^yiy2'  estlnates  of  these  parameters  be  denoted 

—  2  2  —  -  2  2 

by  jf- ,  x„,  a  ,  a  ,  r  and  y, ,  y.- ,  b  ,  s  ,  r 

’  1'  2'  xj'  xg'  x^xg  -^V  yj'  yg  ^1^2 


2. 

”*i 


2  2 


2  2 
t  •  0, 

and  p 

ys  2 

*1*2 

>  and 

p  “  P 

^1 

*2  ye 

against  the  alternative  that  p  ti  n  and  u  A  p  .  Let  the 

the  dtffftre.’ices  bstijean  the  ssuaple  isaas-s  (m  denoted  by  dj^  *  7^  - 
i  «  i.j  9.  fhese  differences,  under  t-hs  cull  hypothesis,  are  cersaU-iy 


dl  mtau  mluss  zero  &s<J  vsriar-cs 


2/1  1  4 

•.‘/7  •  .f  ■ 


V  *  .  A  J.  U 


vters 


/CUM  p.  and  p  csji  b«  estinia'bed  by 


^  (ki-l)8^  (k^-Dv 

8^  »  _  1  *^1  &nd 

Jc^  +  *  2 


(k  -i)a  +  (k  -3 
■*’  i‘'2  ^  ^ 


ki  +  kg  -  2 


1^2  / 

•l*2 


*1*2 


cJt  j5i<*ij  -  ■  '2^ 


tod  «  Is  defined  slnillarly.  To  test  the  hypothesis  that  d, 

7172  3- 

ud  dg  cone  from  populations  with  mans  tero  a^lnst  the  alternative 
that  they  do  not  corns  from;  populations  vlth  means  zero,  ve  fora  the 
expressloo 


^  -2r^ 

*a^dg 


i-r  '» 


<^1 


8 

Sj* 


6:i 


2  2 

In  6.7  a.  and  s^  are  estimates  of  the  variances  of  d,  and 

d^  dg  -1 


dg  and  are  given  by 


T  1*  di  Ktrl  «« 


-/I  1  \  4 

“l^k,  *  '  1-1, 


2  (k^  +  kg  -  2) 

~Z  ^  7/2,  k, +k,,-3 

(k^  +  »g  -  3)  '  i.  2  - 


u,-.'  fJKT.  .a-  .-  ■- 


lAiriol-  sultsbis  ncfJM-lity  asauaiptlcria .  kgsin,  an  teat  la  ra&a* 

Ir-  ths  UBiial  vay  making  u»*  of  the  appropriate  f  tables. 


6.3  aelatlon  hetveen  Trace  B  and  Hotelllag'a 

Returcifjg  nov  to  tiase  B  as  defined  in  Chapter  3  It  t»  « 

g 

exercise  In  algebraic  Kanlpulatlon  to  relate  trace  S  and  T  , 

Recall  tiiat  (xj^,  ^1^'  t  <  1  <  t  are  a  set  of  difference  Tector*. 

One  traoefonnation  taking  (x^,  y^)  into  (a^,  s'>ch  the.t 

k  k  k 

.2,af  .  -  1  and  -  0  is 

3, « j.  1  1  a«j.  i  j,  mi  5  i 


k 


"1 


‘  Ma 


and 


^1  -  T  , 
iSl’'! 


6.6 


i(y  .ilifiilx  i 

Ji^Vi  k  „  *,i,/ 


1=J.  1 


Substituting  6.6  In  the  formula  for  trace  B  and  denoting  £ 

tal 

alatply  by  S  we  obtain 


trace  B  -  +  (Ib^)'^] 


1  ^^l) 


'1  Ek 


DC  > 


E(yi  - 


Exf  ^ 


^  ~  (jhc^}^(EK^y^)''  +  (Exp“(ry^)^  -  I'EXjyjDijIJ'^EXj  +{ lx .J" ( it'' 

EXjI&jEyj  -  (Dc^y^f] 


-  6i>  - 


■ 


Exfily:  -  (Ex,y,)^ 


t«t 


R»en 


»r«ce  B  • 


" 

k5f 

- 

ky 

- 

Ik-l  )s*^  + 

X 

kX^ 

“ 

(k-l)e^ 

+  ky^ 

■ 

(k“l)fl  + 

'  '  xy 

kj^ 

r  ■ 

8  /as 

xy'  X  y 

k®x^((k~l 

)Sy  +  ky^ ] 

+  k' 

yt{k-l)8^ 

+  kX^l  - 

2[(k- 

JSL 


kJtyJk^xy 


xy 


2_2  2  2-2  2  2 

k'^re'  +  kVa^  - 

- _  -  3'S - -2-?*- 


kt(k-l)(s;8^  - 


(k-1) 


1  -  r 


-  +  1 


X  y  8 


froo  foncula  6.5  >»«  see  that 


[k  -  2kr  -iil-  *  k  4 
,  k'  iT  6  s  2 

1  -r  a  X  y  s 

X  ^  y 


Thua 


wane?;  o 


(k-l)-2  1  (k-i)  +  TT 


6v4  Distribution  uf  trace  B  UncLor  ?<ora  a  llty  AoBUsiptlonB 

Tiie  vilstrlbutloo  of  T  Uuuor  dult&uXe  nonaaiiiy  aB&U2iptions 
fur  p  dimer.s!o/ij  and  k  replicatiane  Is  given  by 


K-AAbme  '• 


1 


fCr'idT*' 


B(^- 


t  -1  -• 

k*"irr;  \kT; 

"#  rV/ 


I^t  tr*ce  B  vben  the  urxSijrlylDg  dlstrlbtitlcsi  Ib  nomai  be  denoted 
by  Wjj,  Sine® 

W  -  — i - 5 

"  (h-l)  4^  'T 

for  the  ease  of  two  dlBenslonB,  one  arrives  at 

2 


2  ^ 


dH. 


i  -  W, 


B(V,  1) 


6.9 


-or  the  dl  stributl-/'!  of 


6.5  An  IgilrlcBLl  Sajipllng  ixperlaent 

Per  the  purpose  of  investigating  entplrically  tns  effect  of  non- 
2 

normality  on  Sotelling's  S  ,  a  san^llng  >xperlicont  using  random 
nugibers  and  random  normal  deviates  was  performed.  Forty -eight  sajnples 

of  size  3  were  diawn  from  both  a  two-dimensional  uniform  distribution 
on  ( -i  j  and  from  a  bivariate  normal  distribution  with  zero  mean, 
unit  variance  and  zero  covariance. 

The  sample  size  8  was  chosen  for  several  reaBonc.  It  was  a 
sis.agcablc  size  for  the  necessarj  computations  (Involving  £'  r8ndrji.U m- 
tlons  for  eacii  csjnple)  which  could  be  conveniently  irograniBied  for  the 
electronic  computer  available  (Burroughs  KlOl).  A-so,  It  seems  not 


-  fxS  - 

uijllkely  that  8aii5>les  of  size  8  r  ftfir;  occur  In  actual  problewa  . 
Furthermore,  for  a  aarn}  Is  aa  acali  c;-  Bl7,e  8  one  Sccomes  sonieirtiat 
concerr*^  ai.out  making  t  ie  neccBear}  normality  ai aunpti  one  for  « 

T®  test. 

Forty-eight  sasgilec  were  chosen  becau^i;  it  seemed  that  a  reason¬ 
ably  large  number  of  samples  -«ould  oe  required  to  make  definite  con  ¬ 
clusions  aboi  t  significance  points;  but  at  the  sane  time  forty-sight 
vas  small  enough  to  keep  tbs  amcunt  of  coc^puting  vlthln  reason. 

Sallies  from  a.  rectar.gular  dlstrlbutloo  were  chosen  because  of  the 
coDveclence  In  using  randoai  DUBOers  and  the  desire  t''  Investigate  a 
distribution  co.-isldorabiy  different  free;  the  normal  distribution , 
3a*5ilea  from  a  nonsal  dlstiioutlon  were  chosen  to  both  verify  t)ie 
results  arrived  at  by  direct  mathematical  seans  wd  at  the  same  time 
to  C.1  e-ariy  show  the  contrast  between  frequency  distributions  of  sig¬ 
nificance  points  obtained  from  saaqjles  from  a  non-nomal  distribution 
and  those  obtained  from  saagiles  from  a  normal  distribution. 

An  individual  pair  of  elesmnts  from  any  sample  of  slie  8  was 
considered  as  being  the  observed  difference  of  the  meeaurececti  within 
a  block  from  an  experiment  consisting  of  eight  blocks,  each  block 
coitainlng  the  rerults  of  tne  quantitative  neaBurement  of  two  responses 
from  two  different  trestioents.  Those  differences  (let  thorn  be  denoted 
by  X.  .  V.  )  are  then  renreieh  in  *ccor±tncs  vith  the  trshsfoiiiiiiuu 

i  •  •  1  ■ 

given  by  formul-.  6.0  with  k  «  8.  For  each  6aag)le 

®  2  ®  2 
0  trace  B  -  +  ^iSl^^l) 

was  coapiited  as  wcie  ail  the  127  other  traces  obtained  by  not  changing 


This  rise 


B.nd  changing  the  sign  of  each  pair  of  (a^,  ) . 

a  sat  of  128  quantities  ranging  frees  0  to  5 


with  ao  average  vsdus 


of  a. 


With  the  128  valaee  of  3  tsinoa  B  for  each  sftagsl*,  ob«  can 
piclt  out  the  largest  VBJ.uae  and  then  coaqsare  the  aiapirlcaX  distrlbutloD 
vith  that  expected  under  nornaJLltjr  asaug^Ttiuda.  A  freqiMncy  distribu¬ 
tion  of  ttse  nuidber  of  values  of  W  •>  8  trace  B  which  exceed  the 
point  of  Wjj  oaa  be  coRi.tructad.  Bie  frequency  olstributloee  for 
both  the  sajaplee  fross  a  normal  dlntributloo  and  from  a  rectangular 
dl»trlfcutl<M3  for  the  5.0,  5?. 5,  1.0  and  0.5)6  point*  are  given  in 
Table  6.1.  TS.ble  6.2  give*  tha  actual  pereentega  point*  of  the 
eaplrJ.aal  distributions  obtained  by  coagsutlng  the  mesn  of  tb.?.  distri¬ 
butions  in  Table  6.1  and  converting  to  percentages.  Also  95)6 
coefidenoe  Interval*  are  given  in  Table  6.2.  It  follows  fi^aa  Table 
6.2  that  an  apisroxlnate  rule  for  ccaqsutlng  the  actual  ■l^lficance 
level  sdien  a  test  1*  carried  out  at  the  (a  >  o)  level  on  data 

which  deviate  from  noiiml.tty  to  the  saiae  way  that  the  unlfoiet  distri¬ 
bution  devlatea  from  normality  Is  Blagily 

actual  slgoificance  level  -  l.i  a  e  0.5 
This  lormula  gives  actual  algniflcance  levels  of  1.05,  1.6C,  3*25 
and  6.00  In  place  of  the  observed  O.96,  1.68,  3.24  and  6.02. 


6.6  Adjustment  of  Parameters  in  Slstrihiit.iae  of  Trdcc  * 

In  the  foregoing  example,  one  Is  left  with  the  possibility  of 
altering  the  pareaeterfi  in  the  distribution  of  trace  B  so  a*  to 
Hake  a  more  nearly  exact  significance  test.  This  (»n  be  draie  in 
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6.1 

/requen'jy  dletrlbutlcn  of  upper  a$  tail  baadd  on  lioruiai  theory  for 
*♦6  sauKpiae  from  the  rectanRo^ar  distribution  and  the  nonasl 

dlatrtbutlon 


Ko.  exceeding 
cti  point 

^orsa.1  Sasiples 

. 

Rectsuigui.&r  Samples 

5Jl  2.%  l.o?t  0.5)1 

|^|||||9k|j|||||^^ 

0 

1  14  22 

128 

1 

4  13  17 

1  1  13  2i 

2 

3  11  12  £ 

0  2  16  13 

3 

1  7  8 

0  10  10  3 

W 

491 

1  14  7 

5 

7  12 

4  12 

6 

9  4 

4  6 

7 

6 

12 

8 

7 

9 

9 

6 

9 

10 

3 

4 

11 

2 

4 

Table  6.2 


A'^ual  si^hilflcance  level  of  test  for  saogples  froa  noraal  and 
recta’jRuIar  distribution 


lornsLl  Saagiles 

Rectangular  Saatplet 

lOKinal  $ 

5.0 

li.5 

1.0 

0.5 

5.0 

2.5 

1  A 

aL  S  V 

0.5 

Actual  jt 

5,18 

2.72 

1.05 

0.60 

6.02 

3.24 

1.68 

0.96 

95^  Confideiice 
Interval 

— ^ — - — 

kM 

to 

5.70 

a. 36 

to 

3.08 

0.79 

to 

1.31 

0.4o 

to 

0.80 

— 

5.58 

to 

6.46 

2.94 

to 

3.54 

1.42 

to 

1*9^ 

0.78 

to 

I.l4 

-  VJ 


r.  nuiDb'^r  of  ay  cairinji  c<!jr;t*ln  R:53uir;ptlonc  tj.-i  thi  ;'onii  of  the 

betj»  cLlstj/ibutior,  to  be  fitted.  'i'he  west  generiti  modiei  that  ooa 
can  fit  laakaa  use  of  tiic-  first  r'oar  .ruaidarjts  of  ti-»ce  B.  This 
r«<suire»  t.hc  fitting  of  tha  four  paxarsetars  In  ths  dlstrlbiitlor)  glTCr; 
by 


„  H-1  „  v-l 

f(W)<W  -  C(j  ••  1)  (i  -  g)  »  <  W  <  b  6.10 

vnitra  C  1«  a  constant  r  ch  that 


b^ 

j  f(V)(JW  -  1 
a 


Since  we  linov  that  trace  B  Isa  randoa  variable  betvaen  0  and  1 
It-  «aem«  that  «.  oonveaient  thing  to  do  Is  to  fit  n  and  t  by  the 
first  two  cuBulants.  Itae  problem  then  beccoas  one  of  fi.'idlng  the 
afi  po.lnts  of  a  distribution  of  ths  fom 


du 


6.11 


In  the  case  of  sas^leo  from  a  ncrw-1  dlst.  ibutloo  p:  »  1  and  v  ■  3; 

so  ooe  is  then  required  to  find  the  a  i  points  of  6.11  In  ths 

neighbourhood  of  p  ■  1  and  v  ..  3.  Since  the  irablee  of  the  Inocss- 

plete  Beta  rxictlons  give  only  values  of  p  and  v  by  steps  of  O.J 

ana  since  Interpolation  for  such  small  values  of  p  and  v  la  very 

unsatisfactory,  one  must  approxlaate  these  percentage  points.  An 

“pprcjilmstc  t^cthod  of  fluding  the  necessary  percentage  points  is  given 

In  the  Appendix  togu' her  with  a  chc.rt  (rigure  A.l)  showing  the  i^iper 

0 

5.0,  c’.5,  1.0  and  0,54  curves  of  fiw  over  a  range  of 

p^CL  '  p  p 


Met*  ttjat  Vn'.ausF  t.h-j  i.it-li'i  of  ■i;'»  1 1'9  j.-ftrautatlono  or"  5  is 

l/'i  fo'"  *',’*ry  aempie,  ti'-c"  ratio  .>r  u  to  o  1  is  alvays  i/j.  By 
Bailing  u«e  of  ?igurti  A. I,  or."  o»n  fi.-ni  psi  cc-:it«i«e  polsitfl  to 

apply  to  till"  riiiiilcMi'.t  of  thi?  eeiiaplee  frc,it.s  the  r^t'taiigal&r  dls- 

tributton  .  The  fre^tieflcy  dt strlhutions  ere  glvan  in  Tesle  6..j. 


TVbie  6.3 

Frequency  distributions  of  upyicr  cef,  t«l.l  based  on  adjusting  the 
degree#  of  freedom  In  the  uncooditlofiai  distribution  of  V  by  iieaui# 
of  the  first  two  cuBUiRnts  of  the  condltlonai  dlfitrtbutloo  of  W  for 
esraplea  fros  the  rc ftangui.*Jr  distribution. 


No. 

exceeding 
:rf  point 

Boctangular  Sac,.,  's 

%  2.%  l.(rf  0.% 

0 

.1  k  it 

1 

1  21  29 

2 

1*  16  5 

3 

16  5 

2  19 

5 

3  2 

6 

12  3 

7 

15 

6 

7 

9 

•3 

10 

J 

1 

Table  D.ll 

Actual  signlf ; ennee  level  of  erf  test  for  seaples  from  the  ractangular 
distribution,  adjusting  u  ajid  v  ly  the  flrot  two  cumulsrits 


Nominal  $ 

Actual  $ 

95^  Confidenoff’  lutr^r'ifaJ. 

P.C 

5-16 

k.HR  t.n  5.46 

'.5 

■’.71 

2.V6  to  2.98 

-1  -0 

1.17 

0*^#  to  1  #3* 

0.^ 

O.t'o 

0.51  to  0.75  j 

Ttiivl,!  ''.4  .'.howa  the  filgnifi caiir.c  levei:)  aft-;?  liVJikjr:^:  tlie  tti!,;ust!;k;r;t 

ajid  on"  '«.«  3ee  that  this  ojetho-i  of  ndjiist.'ne-nt  works,  qi-Ute  veil  in 

Ihle  inetor^ee  .  Ti^urc  6-1  ohovs  the  liiUiv  iuu/el  obtained  by 

Belecting  a  value  frctn  the  peiTsutatliari  sUstribution  In  each  of  the 

96  sainpies  which  is  such  tb»t  5*4  of  the  va.l'J!eo  of  6  trace  B  are 

larger  tlicui  the  value  eelected.  The  fact  tliat  the  points  obtained 

frou  the  saaiplea  from  the  rectangular  dietrlbution  wlJen  plotted  against 
8  0  2  2 

E, (a  b  )  (a  function  of  the  variance  of  trace  B  only)  have  a 
p  *•!  *  P  P 

Bleiilar  pettens  to  thoee  obtained  frat  the  saiiqjles  frcsn  the  norjual 
oiatrlbutlon,  Is  Airtiher  Juattflcatlou  for  fitting  the  upper  tali  of 
theoe  eipplrlcal  permutation  dlatributions  obtained  from  the  sample® 
from  the  rectangular  distributions  by  the  first  and  second  cumulanto 

8  o  2 

oTily.  In  fact  tl<e  regression  lines  fitted  through  t  '’p)"'  " 

Sind  8  trace  B  «  5-05»  the  point  corresponding  to  p  «  i  and  v  «  1, 

are  almost  identical  for  the  two  enplrlcal  permutation  distributions. 

A  plot  of  the  deviatiooB  of  the  actual  ^  points  from  those 

based  on  the  fitted  beta  distribution  against  the  sixth  degree  invsjriaat 

k  2  13 

(  E  (a  +  b  y)  shows  that  little  or  nothing  would  be  gained  by  fitting 

P«!l  p  p 

the  beta  distribution  by  the  third  cumulant  also. 

Other  ways  of  adjusting  the  beta  distribution  to  fit  the  upper 
tali  of  the  empirical  distributions  were  also  tried.  However,  these 
siethods  proved  to  be  considerably  inferior  to  the  orie  Just  outlined. 

They  are  presented  here  to  Indicate  ways  of  Improving  the  sigiilflccnce 
levoi  of  Lht  test  when  one  does  not  have  the  tinie  or  inciinetlon  to 
approxlflatc  the  percentage  points  of  the  incomplete  beta  function. 


l>i)rs>*a»dS  UiKjb.  . 


I  ; »!  6  %T.  ft.  1 1  ( 


H  a  fi! 

*  i  '  jt  ' 
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Tb.L>  ‘  1 


ii  V  ;t  i. 


,1  .  str  SS. 


L  f  ii 


^  I  ^ 

yi  y;. 


U  t  t.iv:- 


oppo.r.j4> 


'i;  'TiJSi  tP 


I’iiDjf  A  ^  itra-.B  B' 


t,r«i,'»  B, 


Tflii  IP*  A  ^  0**  t 


Twup  A  1 


•  V  i  a  ' 
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y'-irt.hf r  Methods  for  Afl,)iietlng  SXgnlflcMioe  Levelg 
The  criterion 

*  -  tr«c5  B  «  a  )*■  +  ( 

1  P  1  P 

be  considered  »s  being  distributed  as 

ii 

f(W)dW  -  c(l  -  j)  an  0  <  W  <  a  6.12 

There  are  B»oy  reasoos  for  rejecting  this  distribution  froa  0£W9ldera- 
tioo.  The  main  reason  for  rejecting  Ig  probably  bscsiuse  we  know  that 
the  upper  limit  of  the  range  of  trace  B  Is  1  but  chlg  distribution 
given  by  6.12  has  an  upper  limit  of  a.  Thus  one  might  arrive  at 
a  ctji  significance  point  for  trace  B  in  a  particular  caae  which  is 
outside  of  the  range  (o,  1).  However,  this  distribution  doeu  allow 
one  to  sdter  the  significance  points  with  a  very  small  amount  of  cou- 
putatloc.  To  fit  this  distribution  one  needs  to  find  a,  v  and  c 
such  that 


o  v~l  , 

J  ‘  M(1-  -  -1  ”  17 

0 

a 


f  0  '/{I  -  -iW  »  i  ^ 


A  O  ^  X  I 

-  Ha-  *  ^ 

32  X'  P  P'  '  lb 


0 


ilicoe  e(iiiatioiiS  when  solved  give 


1,  2  ^2,2 


i:(a^  + 

1  ?  i- 

'i'hefce:  i'orKiulas  applied  to  the  *rt5  Baciiilee  frorii  the  w:Jfor.'!i  dlatr;- 
butiOD  give  rise  to  a  set  o!  v'c  ajid  a'a  each  of  whien  are  used 
to  obtain  a  vcSae  of  by  the  fojuiula 

w 

(1  -  _a)  -  a 

a 

or  -  a(i  - 

Orse  then  records  the  nijiritcr  of  values  in  the  randord zatlon  distribu¬ 
tion  of  trace  B  wblcU  siceads  this  The  frequency  distribu¬ 

te  ons  arc  given  In  Table  6.^  and  associated  adjusted  percentage  points 
and  conflfience  Intervals  in  Table  6.6.  It  is  interesting  to  note 
that  the  nominal  5>  and  S,%  slgsilflcitnce  levels  for  the  case  of 
the  Munples  from  the  uniform  distribution  give  rise  to  actual  observed 
significance  levels  of  ^.17%  and  as  opposed  to  the  6.0li 

and  arrl'red  at  witnout  any  adjustment.  However,  the  1^  and 

,%  levels  do  not  have  this  rather  good  igjeement. 

An  alternative  method  Is  to  fix  the  upper  end  of  the  distribution 
at  1  and  allow  the  lower  end  to  v-ndsr.  TIiIe  agiln  i-equlree  the 
use  of  only  the  first  two  eumutsmts  of  B  and  without  further  details, 
the  results  Hre  giver,  ir.  Tablos  6.7  ana  6.6, 

Neither  of  these  two  approaches  gives  rise  to  very  sat! s ' actory 
reii'jl  la .  However  the  first  ol  the  two  does  give  reaaon.able  .i-ejulti; 
at  the  'jja  and  2. fit  levels  at  the  -cost  of  bad  r-esuits  at  the  1$ 
and  0, >,i. levels .  if  ooe  were  interested  In  a  test  .at  the  ‘yfa  ie.'el 


ami  wished  to  rttaXe  an  iy.'pru‘Ve.nw?nt  in  the  accuracy  of  the  actua.1 


Frequency  ui strlbiitions  of  upper  Ofi  t*!).  "pased  or.  »,llf:wSr.g  the 
le.uer  end  cf  tia  ilstiributi™'  of  tr^r*  is  to  be  fixed  at  zer.  .1 
ail,1U3tlng  the  upper  end  ftiid  ths  penaraeter  v  for  Baiardss  from  the 
rectangular  distribution 


Ho.  exceeding 
acft  point 

iftcctangulu  Sax^Xeo 

at  .  1.0^. . 0.?i 

0 

lit  .13 

1 

2  29  9 

2 

13  It  1 

3 

lit  1 

k 

2  lit 

5 

9  5 

6 

U, 

7 

Ik 

8 

8 

9 

3 

10 

1 

Table  6,6 

Actual  olgnlflcanoe  level  of  test  for  saiqiles  from  the  recteui- 
gtlar  dlatrlbutlon,  adjusting  v  and  upper  end  of  dlotrlbutlon  of 

trace  B 


Homlnal  ^ 

Actual  i 

9%  Confidence  Interval 

5 

5.17 

4.67  to  5.47 

2.5 

2.46 

2.22  to  2.70 

1.0 

0.65 

0.51  to  0.79 

0.5 

C  .lO 

0.00  to  0.20 

■'■si’oie  6.7 


uptier  end 

adjusting 


aiBt'ioutioKS  of  upp'ir  njt  tail  b«;!ei5  on  (tllo'/ing  the 
of  tho  dititrlbutlou  of  trace  B  to  b.',  fixed  nt  one  aj)d 
the  luoer  end  and  the  pararooter  v  lor  soiiiplen  from  the 
rectangtslai-  distribution 


Ho.  exceeding 
ert  point 

|||||||||||H||||^ 

0 

1 

1 

1  20  27 

?  17  8 

■> 

lit  9  1 

it 

2  2l* 

5 

9  1 

6 

8  4 

7 

16 

8 

7 

9 

7 

Table  6.8 


Actual  elgnlflcacce  level  of  a$  tent  for  samples  from  the  rectan 
gular  distribution,  adjusting  v  and  lover  end  of  dlstrlbutxon  of 

trace  B 


Moailnal  ^ 

Actual  it 

9%  Confidence  Interral 

5.0 

9.3? 

5.01  to  5.63 

2.83 

2.57  to  3.09 

1.32 

l.l4  to  1.50 

0.75 

0.59  to  0.91 
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eignlJlcance  .'ioveL,  then  thlc  aiethcd  does  offer  a  fairly  Bimple  and 
yrct*tly  rsllativ  iuiBVor.  Little  enn  be  B«id  in  imvor  or  uelng  tiie 
second  of  the  two  siethcds  except  tluat  the  imther  poor  layroyament  In 
the  actual  elgnlfloanoe  level  1»  uaifons  for  all  four  levels  investigated. 

Any  other  method  of  fitting  using  the  beta  distribution  requires 
the  use  of  leore  extensive  tables  than  the  ones  now  available  for  obtain¬ 
ing  the  percentage  points  of  the  beta  dlstrihutions .  If  these  tables 
vers  available  over  a  floe  grid  for  snail  numbers  of  degrees  of  freedom, 
then  one  could  test  the  adequacy  of  fitting  by  only  the  first  two 
cuBulants  versus  the  use  of  all  the  four  'Aisiulants  which  have  been 
tbeoretlcskUy  detemlned  and  can  be  coiq>uted  rather  easily,  Xables 
of  the  necessary  coiiq>utatloos  for  obtaining  the  first  four  cwsulants 
of  trace  B  for  all  aaaq>les  are  given  in  the  Appendix  (lables  A, 3 
and  A. 4). 

6.8  Sugaestlons  for  Further  Investigations 

Although  certain  conclusions  seem  evident  from  the  llsdted  ea^lricnl 
study  done  here,  it  seems  that  In  order  to  fully  Investigate  the  effect 
of  noo>noraallty  on  Botelllng’s  one  needs  to  Investigate  other 

distributions  than  the  rectangular.  A  dlatrlbution  which  would  seem 
to  be  vo.rtb  considering  Is  the  double  ei^ooentlal.  This  would  give 
some  inronnation  on  a  distribution  which  Uffers  from  the  normal  dis¬ 
tribution  in  the  oppostle  way  to  that  Ir  which  the  rectangular  differs, 
fUearly,  there  is  no  end  to  the  number  of  sasqjle  sites  that  one 
mlgtit  use  if;  the  permutation  approach.  However,  results  based  on 
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sasrpies  of  b1i3  6  only  are  certainly  Iniftsiuato  In  diswlng  i;on- 
cluulODB  about  saqblea  of  alzas  other  than  t  without  flrat  In- 
veitlgatlng  what  happens  for  aevaral  dlffo.nsnt  eample  slzeB.  To 
Invi'sstlgate  8asg>le(  much  larger  than  stxe  6  one  neebi  at  least 
a  and.!. lira  site  electroalc  coo^.iter. 

'Hia  proble*  of  Investigating  Hotelling's  1^  wliei'e  thexe  are 
more  than  two  rc  sponsea  being  snasured  uaa  not  treated  here  Itit  a 
snail  sjscwnt  of  uodlft cation  In  the  snthcds  used  abculd  product  the 
necessary  rasolta.  Coi^utatioosaiy.,  aore  than  two  rosponses  should 
create  no  n«v  prdihle.iM. 
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AKPflDIX 

A .  1  Method  or  Approxlfnatlng  the  i-'orciintagc  PoUito  _C5f  the  Beta 
Pletrlbutloii  Over  e.  iialtad  rarige 

W«  require  for  a  «  0.05,  0.085,  0.0.1,  C.OOJ  where 


I 

B(!^jp  v) 


J  (1  -  w)'"^ 


dW 


over  a  reoge  of  v)  such  that  v  ••  lies  between  0.5 

and  S.O.  If  ve  let  1  -  U  -  t  then  A.X  becomes 


X-Vq 


B(ti,  v) 


t*"'  at  -  3 


Kx^itoding  (1  -  t)***^  by  e  Thylor's  series  we  obtain 

— i—  /  ll  -  (M-l)t  +  (>.-l)(ii^)t2  .  (^.i)(l^)(K^)t33tV-]  .  a  ^,2 

»(»*,  v)  ^ 

ttoplaclrg  d  -  fey  t^  and  Integrating  tern  by  tera  A. 2  is  replaced  by 


.v+1 


^V+2 


►  v+3 


~  -  (ii-i)  -^Y+  "<  oBlk*,  v) 

Applying  MevtoD's  nsthcd  for  finding  the  root  of  a  poXynoBlaX, 
can  be  quickly  found  over  the  necessary  range  of  v  ~  3m.  The  accureu:y 
of  tbla  approximation  vas  checked  at  each  end  of  the  range  of  ^  re-, 
fuired  (0.5  end  2.0)  against  the  values  obtained  by  li.neor  Interpo¬ 
lation  In  the  Tables  of  the  Incomplete  Beta  fvsnctioti.  The  approxiise. . 

t4  .M-s  .H  4  -^4  ^4  A4  nf  4^  ^  <k  1  ^  .Met,,,..  W 

bv>  CAW  xticao  w  a  wavuT  e  j,  4  A.w.rwAi «w  va4.|^  wp  a  waa  caA.A.  a  ua  vac 

A  chart,  cor.stitictod  from  the  cor..,)utod  poir;ts  for  p  »  .7,  .8,  .9,  1  .1, 
1.2,  1.3,  i.l,  i.5  arid  2.0  Is  given  Bta  Figure  A.i* 
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K.?  Ofmxiit.  oi.  Tuvlea  A. I,  A.li 

Ji  Boine  or  the  earnpies,  mare  than  kh«  tor;  iargaat  valuea  of  the 
periciitatlon  ul5t.ributlon  era  gtron,  Thar-e  extra  yaiuea  a«;  git'OB 
tecausa  those  saaipleo  poaEsaeed  more  tiian  tea  pcrtiMtotlona  exceeding 
the  ^  normal  theory  elgnificance  point. 
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Tne  V»ii*es  of  th«  Ptnuutatio,!  Dist-rtbutior.  of  (rrsc-is  B) 

h.^.  QAnmiM*  fl  ;v. -5  w*  f  i 


for 

48  of 

8 

froH'  th-  ' 

LJclfers  2i 

strJcatios 

SAKPtl  HO: 

1 

2 

3 

k 

■5 

6 

1 

B 

7.2itO 

5.119 

6-503 

6.961 

6.682 

7.180 

7.286 

T-579 

6.T63 

k.99S 

6.067 

6.866 

6.5£5 

5.713 

6.761 

T.221 

6.^9 

5-990 

9.%1»» 

6-661 

6.525 

5.589 

6.315 

6.020 

6.1»U 

*.9>n 

5.360 

6.3ca 

6.173 

5.5«1 

6.032 

5-888 

V.836 

5-308 

5.710 

5.789 

5.513 

5.670 

5.I37 

6.153 

U..672 

5-158 

5 -116 

5.363 

5.352 

5.362 

5.110 

5.^32 

4.620 

5.112 

5.395 

5. 318 

5.152 

5.289 

9.2S6 

5.357 

k.6n 

5.091 

1-901 

5.259 

1.896 

5.005 

5.271 

5.0*1 

H.itSS 

5.067 

1.733 

5.222 

1.892 

1.979 

fc.992 

5.031  i».*»75 

OAMPLE  NO: 

5-036 

1.638 

9.211 

I.63I 

1.971 

1.987 

9 

10 

11 

12 

13 

11 

15 

16 

7.555 

5.761 

6.305 

6.960 

6-991 

7.277 

7.280 

7.205 

6.130 

5.589 

5-768 

6.137 

6.618 

6.911 

6.955 

6.677 

6.051 

5.535 

5.689 

5.708 

6.512 

6.073 

6.I80 

6.0II 

5.730 

5.516 

5.697 

5-502 

6.071 

6.071 

5.658 

5.866 

5.503 

5.3I8 

5.631 

5.119 

5.719 

5.661 

5.360 

5.622 

5.3I5 

5.321 

1.901 

5.262 

5.670 

5.612 

5.213 

5.385 

5.331 

1.829 

1.895 

5.169 

5.602 

9.618 

5.168 

5.231 

5.153 

1.766 

1.773 

5.106 

5.I91 

5.538 

5.130 

5.219 

1,832 

1.675 

1.^1 

5.091 

5.375 

5.006 

5.108 

5.161 

1.815 

1.385 

1.581 

5.011 

5.095 

5.003 

5.0I0 

5.059 

SAMnLX 

HO: 

17 

16 

19 

20 

21 

22 

23 

2l 

6.563 

6.993 

6.816 

6.612 

6.693 

7.620 

7.172 

6.659 

6.185 

6.802 

6.650 

6.365 

6.321 

6.701 

5.792 

6.212 

6.053 

6.135 

6.021 

6.273 

6 .271 

5-719 

6.172 

5.997 

6. 381 

5-920 

5.861 

6.0^ 

6.k)2 

5.^9 

5.551 

5.81v3 

5.527 

5.513 

5.175 

6.017 

5.10I 

5.271 

5.122 

5.59a 

5.392 

5-302 

.'.Its 

5-712 

5-012 

5.198 

1.971 

5.589 

5.158 

5.158 

5  - 108 

5.5II 

1.889 

5.165 

1.526 

5.I69 

5.087 

5.131 

5.266 

5.013 

1.005 

5.151 

1 .650 

5. I'D 

1.702 

5.10I 

5.062 

1.6il 

1.626 

1.782 

1.396 

1.903 

1,663 

5 .050 

5.019 

1.598 

I.!?*" 

1.396 

T»bi«  A-l  (coijtd.) 


oAHPLF  NO: 


pK 

J-3 

2? 

26 

29 

:;0 

32 

6.723 

b.J'Tii 

7.444 

7.166 

6.607 

7-296 

ti.6ov 

6.681 

fc.364 

6.09^^ 

6.450 

b  ,6^4 

6.221 

6.962 

6.242 

6 . 6.26 

5.S«V 

5-755 

6.341 

6.437 

5  -914 

9 . 520 

6.130 

6.16? 

5.756 

5.5itC 

6.018 

5.519 

5.706 

6.517 

5.725 

6.150 

5.136 

5.527 

5-868 

5.386 

5.794 

5 -ST? 

5.425 

5-754 

5.039 

5.368 

5.8'-' 

5.087 

5-340 

5-749 

5.424 

^ ,  6b6 

5.03^1 

5.166 

5.579 

^*.953 

5.296 

>.;87 

5.062 

i;  ■  5(0 

5.023 

5.0i*5 

4.904 

4.879 

5.145 

5.395 

4.901 

5.433 

Six; 

it.8a8 

4.5S3 

4.871 

4  .Q44 

5. 

4.884 

5.209 

it. 786 

4.602 

4.573 

4-799 

4.744 

5.209 

4 .876 

5.201 

5.184 

5.171 

SAMPLE 

SO: 

33 

34 

35 

36 

37 

33 

39 

4o 

7.118 

7.807 

6.486 

6.682 

7. '>24 

7-251 

7 .262 

7.655 

5.985 

6.346 

6.375 

5-925 

6.601 

6.627 

6.598 

6.322 

5.978 

6.067 

6.328 

5-790 

6,282 

6.153 

5.669 

5.960 

5.965 

5.879 

6.095 

5.719 

6.178 

5.937 

5-470 

5.»j»6 

5.929 

5.737 

6.078 

5.6J5 

5.885 

5-690 

5.247 

5-732 

5.859 

5.449 

5-964 

5.631 

5.139 

5.493 

5-1.79 

5.728 

5.306 

5.4oo 

5.509 

5-236 

5.142 

5-282 

5.011 

5.486 

4.809 

5.U7 

5.287 

4.546 

5.015 

5.235 

4.969 

5.389 

4,625 

4.887 

5.051 

4.531 

4.946 

5.189 

4.967 

5.369 

4.348 

4.731 

4.930 

4.462 

4-772 

4,787 

4.720 

5.206 

5.153 

SAHPUS 

SO: 

41 

42 

43 

44 

45 

46 

47 

48 

7.553 

7.491 

7.142 

6.786 

7.176 

6.751 

6.448 

7.179 

7.016 

6.355 

6.909 

6.435 

7.rx)2 

6.107 

6.293 

6.351 

6.972 

5.662 

6.569 

5.935 

6.416 

6.064 

6.165 

6.151 

6.354 

5-694 

5.630 

5.616 

6.294 

5-704 

K.  r 

6.227 

5.930 

5.713 

5.615 

5.6^ 

5.580 

5-202 

W 

5-384 

^.V>(U 

5.511 

5.059 

5.008 

?.3S6 

5.162 

5.266 

5-047 

5.444 

5.574 

5.061 

5.344 

5.461 

4 .948 

5.056 

4.842 

5.4ii 

5.538 

4.805 

5.263 

5.317 

4.720 

4.919 

4.668 

5.233 

5.3&I 

4.736 

5-247 

5.268 

4.645 

4.887 

4.652 

4.900 

5 . 346 
5.06'. 

4.698 

5.182 

4.615 

4.618 

4.884 

4.425 

-  82  - 


Tiibl*  A. a 


The  lArcest  Value#  of  tbe  PerButatlon  Dl»t  rihutloa  -  f  &(t.-act  B) 
for  Wi  afsgjle#  of  alM  8  fro*  tbr?  Bcrsaal  DIstriButtoe 


SAMPUC 

so.' 

> 

2 

3 

k 

5 

6 

7 

8 

C.I65 

5.885 

7.323 

5.983 

7.034 

6.939 

6.6?7 

7.245 

:i.727 

5.172 

6.811 

5.729 

6.559 

6.858 

6.731 

6.144 

5.179 

^.7-9 

6.0^9 

5.653 

6.363 

6.8*51 

6.530 

5.84o 

5.115 

k.Wi2 

5.780 

5.467 

6.072 

6.444 

5.894 

5.730 

5.0l^3 

a.iiW 

5.763 

5.063 

6.036 

6.148 

5.770 

5.685 

5.001 

%.091 

5.515 

4,892 

5.513 

5.527 

5.743 

5.549 

i.737 

3-8£5 

5.^78 

4.785 

5.237 

5.423 

5.458 

5.368 

^.535 

3.857 

5-3i'^ 

4.650 

5-231 

5.129 

5.208 

5.175 

*^.357 

,3.855 

5-287 

4.6?3 

5.157 

4.675 

5.072 

5.166 

‘>.297 

3.aj9 

>*.971 

4.621 

4.889 

4.585 

5.034 

5.102 

3Ai<PLJE  IWi 


9 

10 

11 

12 

13 

14 

15 

16 

7.134 

6.783 

5.726 

6.147 

6.520 

5.889 

6,693 

5.606 

6.678 

6.Tl3 

5.583 

5.738 

6.456 

5. 568 

5.962 

6.043 

6.287 

5.470 

5.730 

6.373 

5.535 

5.851 

5.628 

5.827 

5.318 

5.290 

5.340 

5.978 

5.343 

5.756 

5.410 

5.794 

5.076 

5.182 

5-U6 

5.958 

5.3U 

5.739 

5.241 

5.509 

4.901 

5-114 

5.016 

5.565 

5.266 

5.588 

4.072 

5.399 

4,733 

4.837 

5.015 

5.545 

5.157 

4.976 

4.835 

5.305 

4.781 

4.704 

4.982 

5.289 

5.061 

4.967 

4.732 

5.128 

>*.760 

4.676 

4.862 

4.955 

5.068 

4.676 

4.667 

5.055 

SAMFU:  ! 

4,309 

?0: 

4.558 

4.781 

4.724 

4.776 

4.763 

4.581 

17 

18 

19 

20 

21 

2? 

23 

24 

6.140 

7.035 

6.370 

7.673 

6.776 

6.728 

7.021 

6.047 

5.7^8 

6.9H 

6.130 

6.359 

6.605 

5.727 

5.631 

6. roe 

5^816 

6.215 

5.777 

6.051 

6.228 

5.689 

5.552 

6.267 

5-566 

6-133 

5.698 

5.858 

5.972 

5.190 

5.4^ 

6.o44 

5.499 

5.402 

5.553 

5.825 

5.2^^ 

5.093 

5.416 

5.544 

5.375 

5.115 

5.114 

5.697 

5.225 

5.071 

4.927 

5.499 

4.;?65 

5.043 

5.105 

5.604 

5-122 

4.865 

4  .901 

5.448 

4,92a 

4.979 

5.060 

5.486 

5.049 

4.703 

4.845 

5.2^7 

4.904 

4.9S6 

5.103 

4.d4S 

4.439 

4  J 

5.231 

5.096 

4.7^ 

4.592 

'j  .avt 

4.657 

4 . 5'?'2 

4.306 

-  8‘*!  • 


Tabie  A. 2  (contd.) 


SAMd^JI 

HO: 

?5 

26 

27 

28 

6.633 

7.699 

6.436 

6.459 

6 . 387 

6.481 

6.253 

6.014 

5.212 

6.267 

5.920 

5.613 

5.106 

5.758 

5.595 

5.324 

^.021 

5.316 

5.473 

5.235 

U.9S0 

5.113 

4.833 

5.036 

‘•.978 

4.966 

4.822 

4.896 

4.732 

4.783 

4.690 

4.722 

4.703 

4.583 

4.607 

4  .(^6 

4.630 

4.439 

4.481 

4.598 

SAMPLE 

NO: 

33 

34 

35 

36 

6.161 

7,139 

7.12? 

7.234 

5.465 

6.775 

6.671 

6.463 

5.044 

6.618 

5.980 

6.329 

4.949 

6.030 

5.77a 

6.173 

4,948 

5.911 

5.770 

6.0)^ 

4.883 

5-756 

5  741 

5.588 

4.766 

5.506 

5.198 

5.4o4 

4.520 

5.329 

5.331 

5.CO5 

4.42y 

5.122 

4.987 

4.764 

4.230 

5.093 

4.785 

4.722 

SAMPIS 

NO: 

4l 

42 

43 

44 

7.256 

6.631 

6.838 

6.690 

6.04i 

6.074 

6.737 

5.956 

5.827 

5.520 

5.510 

5.903 

5.740 

5.270 

5.432 

5.773 

5.624 

5.245 

5.125 

5.542 

5.292 

5.091 

5.084 

5.410 

5.122 

5.007 

4.838 

5.302 

5.011 

4.747 

4.816 

5.249 

4.948 

4.573 

4.806 

5.055 

4.917 

4.453 

4.703 

4.928 

29 

30 

31 

32 

6.352 

7.049 

6.11.6 

6.528 

5.6^ 

6.350 

5.800 

6.452 

5,472 

6.227 

5.275 

6.206 

5-428 

6.191 

5.2X4 

6.052 

5.366 

6 .04d 

4 .952 

5.513 

5.316 

5.368 

4.942 

5.229 

5.224 

5.230 

4.010 

5.172 

5.194 

5.182 

4.^ 

5.138 

5-007 

4.941 

4.418 

4.960 

4.963 

4.795 

4.373 

4,919 

37 

38 

39 

4o 

6.936 

7.014 

5.528 

6.639 

6.134 

6.618 

5.313 

6.544 

6.079 

6.6U 

4.958 

6.245 

5.977 

5-623 

4.613 

4.9^ 

5.728 

5.557 

4.749 

4.965 

5.337 

5.396 

4,522 

4.937 

5.252 

5.219 

4.508 

4.914 

5.007 

5.012 

4.498 

4.875 

4.9.34 

4.748 

4.482 

4.580 

4.956 

4.611 

4.475 

4.485 

45 

46 

47 

48 

6.585 

6.028 

7.175 

6.?67 

6.219 

5.885 

6.127 

5.949 

6.o4o 

5.848 

5.745 

5.467 

5.949 

5.803 

5.7D2 

5.284 

5.749 

5.13? 

5.489 

4,675 

5.2!iO 

4.959 

5.486 

4.599 

5.131 

4.931 

5.438 

4.541 

5.002 

4.070 

5.342 

4.489 

4.956 

4.346 

5.182 

'•.303 

4.736 

4.654 

5.073 

4 .242 

-  ah  - 

I'&ble  A.  3 


VeJ.u«s  rf  thi  S.’X  Iuv»ri*iit»  for  Coaiputlng  the  Flrut.  Fsur  CiaBulants 
of  tr&c«f  B  for  Vi  saapl#*  of  «l**  8  froo  tb*  Uftifons  Distribution. 


No 

P  !’ 

£(£,*+bp'*  {E( 

P  p' 

E"(a  »  +b  b  )4 
P  9  P  9 

1 

.5815 

.2000 

.0795 

.3382 

.06.27 

.1251 

E 

1.0525 

.6624 

.4409 

l.’,077 

.0680 

.4912 

3 

.7917 

•  3053 

.S12: 

.6268 

.0705 

.2465 

.6057 

.2015 

.0711 

.3669 

.0753 

.1346 

3 

.5752 

.1804 

.C.60I 

.3308 

.0711 

.1290 

6 

.7481 

.3367 

.1704 

.5597 

.0813 

.2162 

7 

.6282 

.2303 

.0941 

.3947 

.0663 

.1433 

8 

.5608 

.1706 

.0^ 

.3145 

.0893 

.1432 

a 

.6521 

.2413 

.0951 

.4253 

.0935 

.1797 

10 

.9124 

.5164 

.3245 

.8326 

.0521 

.3117 

11 

.8126 

.3907 

.2010 

.6603 

.0689 

.2543 

jl2 

.7157 

.3137 

.1561 

.5122 

.0675 

.1959 

13 

.6046 

.2015 

•0703 

.3655 

.0806 

.1489 

lA 

.6067 

.2054 

.0741 

.3681 

.0864 

.1514 

13 

.6124 

.2068 

■0756 

.3751 

.0908 

1637 

IC 

.6165 

.2074 

.0726 

.3000 

.7oJa 

.0782 

.1455 

17 

•0391 

.45^ 

.2793 

.0732 

.2S)48 

18 

.5092 

.1930 

.0674 

.3472 

.0717 

.1372 

19 

.6540 

,22  V 

.0787 

.4277 

.U68 

.1947 

ao 

.6566 

.2397 

.0933 

.4-13? 

.0818 

.1678 

21 

.7055 

■  3082 

.1487 

.4>-78 

.0648 

.2027 

22 

.6013 

.2144 

.0866 

.3615 

.0661 

.1305 

23 

.8533 

.4637 

.2763 

.7281 

.0662 

.324? 

2k 

.8317 

.4413 

,2632 

.6917 

.0539 

.2453 

25 

.7226 

.3009 

-137? 

.5221 

.1048 

.2373 

26 

.7457 

.3-182 

.1424 

.5561 

.1070 

.2463 

27 

.715^^ 

.3392 

.19U 

.5117 

.0683 

.1625 

26 

.73Tr- 

.3526 

.1967 

-5434 

.0582 

-1935 

29 

.6756 

.2469 

.0931 

.4564 

.0670 

.1793 

30. 

.5354 

.1532 

.0466 

.2866 

.061(9 

.1110 

31 

.7084 

.2929 

.1359 

.5018 

.0855 

.1929 

32 

.5069 

.1772 

•0547 

■  3445 

.0653 

.1400 

33 

.6639 

.2326 

.0840 

.44o8 

.1309 

.2146 

3** 

.6279 

.2272 

.0899 

.3943 

.0958 

.1646 

35 

.6688 

.2670 

.U89 

.4472 

.0736 

.1808 

36 

.7833 

.3418 

.1541 

.6116 

=1543 

if 

.5701 

.1707 

.0524 

.3251 

.1346 

38 

.&:;22 

.2132 

.0767 

.3071 

.o8r2 

.1614 

39 

.3730 

.5330 

■  3?5i6 

.7622 

.0481 

.2737 

lo 

.594'9 

.1945 

.0680 

•  35.16 

.0757 

.1336 

4l 

.52^*9 

.1483 

.0436 

.2808 

.0719 

.1104 

42 

.6270 

.2226 

.064-3 

.3932 

.0779 

.1525 

43 

.6553 

.2606 

.1138 

.4294 

.0921 

.187? 

44 

.6446 

.2277 

-0<358 

■  4155 

.0920 

.1739 

45 

.6092 

.21t8 

.0196 

.3712 

.0643 

.1514 

46 

.7360 

.3058 

.1355 

.5417 

.1108 

.2312 

47 

.6702 

.2387 

.0883 

.4492 

.1038 

.1915 

48 

■  7330 

.3230 

.1»6 

-5373 

.1004 

.2452 

-  05  - 


Table  A.i 


Values  of  the  Six  Invariant*  for  Coaiputing  the  First  Four  Cuwilants 
of  trues  u  for  to  of  6  rr*»m  ti»  nOrcaal  rietri buiion . 


Sanple 

E(a*+b*)^ 

E{a*+b“)3 

ECaW)** 

[E(a“+b»)^J^ 

E"(a  a  4^6  b 

E"(a  b  -b  a  ) 

no 

P  P 

P  P 

P  P 

p  p 

P  9  PI 

P  <1  P  't 

1 

1.001? 

.6292 

.4343 

1.0031 

.0513 

.3892 

2 

1.2644 

.9405 

.7237 

1.5987 

.0462 

.7242 

3 

.6550 

.2643 

.1205 

.4290 

.0719 

.1802 

k 

•  773? 

•  35?9 

.ir96 

.5983 

.0688 

.2346 

5 

.6607 

.2614 

.U35 

.4.365 

.0870 

.2002 

6 

.6001 

.2049 

.0771 

.3601 

.0763 

,l4o4 

7 

.6387 

.2378 

.0^ 

.4079 

.0708 

.1534 

6 

.6303 

.2165 

.0791 

.3973 

.0871 

.1567 

9 

.6341 

.2380 

.0886 

.4021 

.0834 

.1665 

10 

.8724 

.4909 

.3028 

.7611 

.0429 

•2o20 

u 

.8627 

.4267 

.2232 

•7442 

.070-6 

.2704 

12 

.6i4o 

.3884 

.2002 

.662? 

.090.:’ 

.2636 

13 

.6483 

,2293 

.0643 

.4203 

.0^ 

.1541 

Ik 

.7644 

.3229 

.1441 

.5343 

.095*3 

.2350 

15 

.7799 

.3872 

.2166 

.6003 

.o;>57 

.2396 

16 

.9019 

.4973 

.2979 

.8134 

.0638 

•  3311 

17 

.8231 

.4168 

.2267 

.6775 

.0518 

.2739 

18 

.5479 

.1606 

.0495 

.3002 

.0^ 

.U8l 

19 

.9212 

.5494 

.3710 

.8486 

.0548 

.2890 

£0 

.6870 

.3060 

.1596 

.4720 

.0647 

.1725 

21 

.7018 

.2751 

.1132 

.4925 

.0910 

.2035 

22 

.6639 

.2546 

.1051 

.4407 

.0686 

.1671 

23 

.8282 

.4496 

.2842 

.6859 

.0844 

■2357 
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